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Zusammenfassung 

Durch Ankopplung eines Helmholtz-Resonators an ein mechanisches Schwingungsgebilde laBt 
sich unter bestimmten Bedingungen der Kopplung, Abstimmung und Dampfung im Hinblick auf 
die Anwendung bei Telephonen und Mikrophonen ein relativ resonanzfreies UbertragungsmaB 
erzielen. Optimale Bedingungen fiir Frequenzbandbreite erhalt man dann, wenn die Abstrahlung 
oder Aufnahme des Schalls vom Schwingungsgebilde selbst erfolgt, d. h. der Helmholtz-Resonator 
nicht im Schallweg liegt. 


Summary 

By coupling a Helmholtz resonator to a mechanical system in vibration it is possible under 
certain conditions of coupling, tuning and damping to attain a transduction coefficient relatively 
free from resonance, for application to telephones and microphones. Optimum conditions for band 
width are attained when the radiation or pick-up takes place from the mechanical system itself 
without the Helmholtz resonator forming part of the sound path. 


Sommaire 

En couplant un résonateur de Helmholtz 4 un oscillateur mécanique, on peut obtenir, dans 
certaines conditions de couplage, d’accord et d’amortissement, une transmission relativement 
exempte de résonance, en vue des applications aux téléphones et aux microphones. Les conditions 
optima, en ce qui concerne la largeur de la bande de fréquence, sont réalisées si |’émission ou la 
réception du son se font par l’oscillateur lui-méme, c’est-a-dire si le résonateur de Helmholtz ne 


se trouve pas dans le trajet acoustique. 


Einfiihrung 

Bei vielen akustischen Wandlern kommt es 
darauf an, ein sehr breites Frequenzband verzer- 
rungsfrei abzustrahlen oder aufzunehmen. Das 
14Bt sich nur mit Schwingungsgebilden erreichen, 
deren Resonanzstellen unterdriickt sind. In vielen 
Fallen ist die Dampfung eines akustischen Wand- 
lers nur méglich, wenn man ein weiteres, stark 
gedampftes Schwingungssystem ankoppelt. 

Die iibliche Behandlung gekoppelter Schwin- 
gungsgebilde beschrankt sich bevorzugt auf die 
Eingangsimpedanz und auf relativ kleine Damp- 
fungen der Einzelgebilde (z. B. [1]). Solche Be- 
trachtungen sind zwar in vielen Fallen, speziell 
fiir akustische Sender, von grofer Bedeutung, 
doch sind sie nicht erschépfend. 

Es sollen deshalb einige Beispiele behandelt 
werden, die fiir die Dimensionierung von Breit- 
bandwandlern wichtig sind. 


1. Einfaches Schwingungsgebilde 


Ein mechanisches Schwingungsgebilde mit 
einem einzigen Freiheitsgrad (z. B. eine Telephon- 
membran), das durch die GréBen Masse m, Steife s 


und Reibungswiderstand r festgelegt ist, kann 
bekanntlich durch ein elektrisches Ersatzbild 
nach Abb. | dargestellt werden [2]. 


et 

Abb. 1. Elektrisches Ersatzbild eines Schwingungsgebildes 
mit einem Freiheitsgrad (elektrische GroBen durch 
mechanische angegeben). 


Die Schnelle & ist dann 


é=F/W, (1) 
wenn die antreibende Kraft 
F=F,-e (2) 
und der mechanische Scheinwiderstand 
; s 
W—r+j(om—*) (3) 


ist. Die Eigenfrequenz eines solchen Schwingungs- 
gebildes ist durch 

wm=s (4) 
gegeben. Ist ein solches Schwingungsgebilde Glied 
einer Ubertragungsanlage, so mu zur Kenntnis 
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der Ubertragungseigenschaften der Ausschlag 
berechnet werden, der mit der Schnelle € durch 


F=jog (5) 
verkniipft ist. 

Bei der Ubertragung eines breiten Frequenz- 
bandes stért die Amplitudeniiberhéhung in der 
Eigenfrequenz des Schwingungsgebildes. Erst bei 
einem sehr grofen, technisch meistens schwer 
realisierbaren Reibungswiderstand verschwindet 
die Resonanziiberhéhung. Gleichzeitig wird aber 
die obere Grenzfrequenz des zu tibertragenden 
Bandes herabgesetzt. 


2. Gekoppeltes Schwingungssystem 


WENTE und Tuuras [3] haben gezeigt, dab 
wesentlich bessere Ubertragungseigenschaften zu 
erreichen sind, wenn man unter bestimmten Vor- 
aussetzungen an ein relativ ungedimpftes Schwin- 
gungsgebilde ein stark gedampftes ankoppelt. Sie 
haben sich aber dabei auf die Berechnung des 
Ausschlages des einen Schwingungsgebildes (&m) 
beschrankt. Es soll deshalb dieses Schwingungs- 
system einer allgemeinen Betrachtung unter- 
worfen werden. Sein mechanisches Bild ist in 
Abb. 2, das elektrische Ersatzbild, nach HAEHNLE 
[2] Analogie 2. Art, in Abb. 3 gezeigt. Das me- 
chanische Schwingungsgebilde besteht hierbei aus 
Masse mm, Steife sw und Reibungswiderstand rm. 
Zwischen der Masse my und einem unendlich fest 
gedachten Widerlager greift die erregende Kraft 
Fm an. Das angekoppelte Schwingungsgebilde 
besteht aus einer Feder mit der Steife sz, und 
einer Masse mr, und diese ist tiber einen Rei- 
bungswiderstand rz, gegen das unendlich fest ge- 
dachte Widerlager gedampft. 

Im Hinblick auf einige Anwendungen soll die 
Aufgabe die Berechnung folgender GréBen sein: 


a) Ausschlag y der Masse mm und 


b) Ausschlag &, der Masse m, gegen das unend- 
lich feste Widerlager bei konstanter erregender 


Kraft Fy. 


Abb. 2. Gekoppeltes mechanisches Schwingungsgebilde. 


Us AY fs 


Fy Us é, 


Abb, 3. Analoges elektrisches Schaltbild zu Abb. 2. 
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Im elektrischen Analogon (Abb. 3) sind die 
gleichen Bezeichnungen des mechanischen Bildes 
(Abb. 2) beibehalten. Es entspricht also die er- 
regende Kraft Fy einer elektrischen Spannung 
und die Schnelle £ dem durch das ganze System 
flieBenden Strom, die Schnelle &; dem Teilstrom 
durch die Induktivitat (mz) und den Wider- 
stand (rz). 

Eine praktische Ausfiihrung besteht darin, daB 
man an eine Membran eines akustischen Wand- 
lers, z. B. an eine Telephonmembran, einen Helm- 
holtz-Resonator ankoppelt. Die Membran wird 
durch eine elektrische Kraft erregt. Der Helm- 
holtz-Resonator besteht aus einem an der einen 
Seite der Membran angrenzenden Luftvolumen 
V_, mit seiner Steife sy und einer Offnung mit der 
Masse mz und dem Reibungswiderstand rr. 

Die akustischen Widerstande (kg/s m*) des Helmholtz- 
Resonators miissen in mechanische (kg/s) umgerechnet 
werden, damit die Rechenoperationen mit gleichen Dimen- 
sionen der Einzelglieder des Schwingungsgebildes durchge- 
fiihrt werden kénnen. So wird die mechanische Steife sy des 
Helmholtz-Resonatorvolumens Vz 

ce o-S 
n= (6) 


solange die Wellenlange groB zu den linearen Dimensionen 
des Volumens ist. @ ist die Dichte, c die Schallgeschwindig- 
keit der Luft, Syy die effektiv strahlende Flache der an den 
Helmholtz-Resonator angekoppelten Membran. Entspre- 
chend ergibt sich fiir den mechanischen Scheinwiderstand 
der Offnung des Helmholtz-Resonators, z. B. eines engen 
Rohres 


4 
W, =r, + jom,y =2L (s + jo ze), (7) 


wenn die Linge L der rohrférmigen Offnung klein zur 
Wellenlinge und R der Radius des Rohres ist. « ist der 
Reibungskoeffizient der Luft. (Akustische Widerstinde sind 
in [4] auch fiir andere Offnungen kurvenmafig zusammen- 
gestellt.) 

In der Praxis wird haufig statt eines engen Rohres eine 
mit einem Gewebe iiberzogene Offnung benutzt [5], [6]. 
Ein gewiinschter Scheinwiderstand Wy, laBt sich leicht 
durch Variation der GréBe der Offnung und der Webart des 
Stoffes erreichen. 


Mit der Ankopplung eines Helmholtz-Resona- 
tors an eine Membran ergeben sich dann zwei 
Méglichkeiten: 

a) Die Membran kann die Schallenergie direkt 
abgeben oder aufnehmen, d.h. der Helmholtz- 
Resonator liegt nicht im Schallweg (Abb. 4). 
Aufgabe ist also die Berechnung von éu 
bzw. Em. 

b) Die Membran kann die Schallenergie durch 
den Helmholtz-Resonator hindurch abgeben 
oder aufnehmen (Abb. 5). Aufgabe ist hier 


also die Berechnung von &, bzw. &1. 


In beiden Abb. 4 und 5 ist das Schema eines 
Telephons gezeigt, das in einer schallharten 
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Abb. 4. Abb. 5. 


Mit Helmholtz-Resonator 
gekoppelte Membran, 
Schallabgabe von der 

Membran. 


Mit Helmholtz-Resonator 
gekoppelte Membran, 
Schallabgabe von der Off- 
nung des Helmholtz-Reso- 
nators. 
Druckkammer einen Schalldruck p bei Erregung 
der Membran durch eine elektrische Kraft Fy 
erzeugt. Hierbei hat das Volumen des Telephon- 
inneren V7 die Steife sv, das der Druckkammer V, 
die Steife sy, solange die Dimensionen beider Vo- 
lumina klein zur Wellenlange sind. Der jeweilige 
Wert der Steife sr bzw. sy) berechnet sich aus 
Gl. (6). Die vollstandigen elektrischen Ersatz- 
bilder beider Anordnungen (Abb. 4 und 5) sind 
in den Abb. 6 und 7 wiedergegeben. Nun ist aber 
im Hinblick auf die spater zu behandelnden Bei- 
spiele die Steife der Membran sm und die des an- 
gekoppelten Helmholtz-Resonators sz; von glei- 
cher Gré®enordnung. Weiter ist sowohl das Vo- 
lumen des Telephoninneren als auch das der 
Druckkammer gro im Verhdltnis zu dem des 

Helmholtz-Resonators, so daf gilt 


ST, Sy < Si, SM. (8) 

Es gehen also die elektrischen Ersatzbilder der 
Abb. 6 und 7 in das der Abb. 3 iiber. 

Der in der Druckkammer vom Volumen V, 

erzeugte Schalldruck ist bestimmt fiir a) durch 


2 
Pw = = - Su: Em, (9a) 
0 
xe co 2 
fir b) durch py) = Vi Sr an: (9b) 
0 
Ew sy My Thy PM i PR 


Abb. 6. Vollstandiges elektrisches Ersatzschaltbild zu Abb. 4. 


sz My sy Pm m, ry 


"sy Le 
Fuy 1h) she é 
L 


Abb. 7. Vollstandiges elektrisches Ersatzschaltbild zu Abb. 5. 
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St ist die strahlende Flache 
Resonatoréffnung. 

Aus dem elektrischen Ersatzbild der Abb. 3 
1aBt sich die Schnelle £4 bzw. &, ablesen. Die 
Berechnung der Bewegungsamplitude €u bzw. &1 


erfolgt dann iber Gl. (5). 
a) Fir die Schnelle é ergibt sich: 


der Helmholtz- 


: Fu Fu 
Wa Le Wee so 
iy aw, 


worin Wy Gl. (3) entspricht, 
W,.=—j—, W.=rir+jomz, (lla) (11b) 


und schlieBlich W,-, zur Abkiirzung des Nenners 
benutzt wird. 

Der Betrag der Schnelle £m hat unter Vernach- 
lassigung der Widerstande rm und rz, folgende 
Extremwerte: 


Eine Nullstelle bei 
W.+W1=0 oder w?m_p—st=0 (12) (13) 


d. h. dann, wenn die Frequenz gleich der Eigen- 
frequenz des fiir sich schwingenden Helmholtz- 
Resonators ist. 


Zwei Unendlichkeitsstellen bei 
Wu (W.+W1)+W.Wi=0 (14) 
oder (w?mm— sm) (w?mz—sL) —w* mz sz = 0. (15) 


Bei endlichen Reibungswiderstanden ry und rz 
hangt die Ausbildung der Extremwerte von ver- 
schiedenen Parametern ab: bei gegebenen Eigen- 
schaften des einen Schwingungsgebildes (z. B. 
einer Telephonmembran) ist ein optimaler Uber- 
tragungsverlauf |&u/Fm| von der Frequenz zu 
erzielen, wenn man ein bestimmtes Steife- und 
Massenverhaltnis sz/sm und mz/mm und einen 
relativ groBen Reibungswiderstand rz, wahlt. 

WENTE und Tuuras [3] berechneten die Mem- 
branausschlage éy fiir folgende Werte (Einheiten: 


cgs): 


NG lash 2,5 1,25 
Bley 6,25 1,25 0,5 3,45 
Ct 12,5 1,25 0,5 2,95 


Benutzt man als charakteristische Abkiirzungen 
die Eigenfrequenzen Gl. (4) und die Dampfungen 


TO r 
peak. wet 


a 16 
s om 16) 
der fiir sich schwingenden Einzelgebilde mm, sm, 
rm und mz, Sx, rt, so ergeben sich fiir diese drei 


Falle folgende GréBen: 
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Die Gré8en der Einzelelemente fiir die hier durch- 
gerechneten Beispiele wurden in Anlehnung an 
experimentelle Muster gewahlt: 


Die Berechnung dieser Beispiele kann relativ 
schnell durchgefiihrt werden, indem zunachst der 
Leitwert 1/W,-; und damit die auf die Kraft be- 
zogene Schnelle éu/ Fu mit einem logarithmischen 
Rechenschieber fiir komplexe Zahlen nach REIN- 
HARDT [7] bestimmt wird. 

Fir die vier Falle sind die Ortskurven in Abb. 8a 
wiedergegeben. Bei tiefen Frequenzen — als Para- 
meter ist w/wm eingetragen — eilt die Schnelle & 
der erregenden Kraft Fu um 90° voraus, im 
Resonanzgebiet sind beide Gré®en in Phase und 
bei hohen Frequenzen strebt die Schnelle einer 
Phasennacheilung von 90° zu. Aus den Absolut- 
betragen des Leitwertes laBt sich durch Multipli- 
kation mit wm/@ der Absolutbetrag der relativen 
Membranbewegung £y/F'uermitteln (vgl.Abb.8b). 
Von der Frequenz w/wm optimal unabhangige 
Membranausschlage &y/F'm sind also etwa durch 
die Falle (2) und (3) gegeben. Durch die An- 
kopplung des Helmholtz-Resonators wird die 
obere Grenzfrequenz des Ubertragungsbereiches 
um etwa 50% erhéht, wenn man im Vergleich 
dazu ein relatiy stark gedimpftes, sonst aus 
gleichen Elementen sw und my_bestehendes 
Schwingungsgebilde mit dy =1 betrachtet (vgl. 
dazu Kurve 5, Abb. 8b). 

: Gema8 G1.(5) kann man die Werte der Leitwertortskurve 
Ey/F mM mit —jwy/w multiplizieren, um die Ortskurve der 
Membranbewegung &);/ Fy zu erhalten (vgl. Abb. 8c). Hier 
beginnt die Ortskurve mit Phasengleichheit zwischen Be- 
wegung und Kraft. Der Absolutbetrag des jeweiligen Orts- 
vektors iiber der Frequenz aufgetragen stellt natiirlich 
ebenfalls die gesuchte Ubertragungskurve |&y/Fy| dar. 

b) Wird bei dem gleichen Schwingungssystem 
der Schall von der Offnung des Helmholtz-Reso- 
nators abgegeben (vgl. Abb. 5), so hat man die 
Schnelle é, bzw. den Ausschlag &, zu berechnen. 

Aus dem elektrischen Ersatzbild der Abb. 3 
ergibt sich 

Fy 


i 


Woes (1+ W/W)” 


(17) 
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Jm(&y/Fy! 


! 
~ 


Abb. 8. (a) Ortskurve der relativen Membranschnelle 
Eu/F'.- 
(b) Ubertragungskurve |£)4/F' | =f (w/oa). 
(c) Ortskurve der relativen Membranbewegung 
Em/F.- 


Diese Gleichung erweist sich fiir die Berechnung 
besonders zweckmaBig, da bereits die unter a) 
behandelten Werte des Leitwertes 1/W,.. be- 
kannt sind. 

Setzt man W,., von Gl. (10) in Gl. (17) ein, so 
laBt sich der Kurvencharakter fir é,/Fi als 
Funktion der Frequenz erkennen. Hier bildet sich 
keine Nullstelle aus, wahrend die beiden Unend- 
lichkeitsstellen die gleichen wie bei Fall a) sind. 
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Der Rechengang ist dem unter a) geschilderten 
ahnlich. Man erhalt zunachst die relative Schnelle 
é,/Fu an der Offnung des Helmholtz-Resonators 
(vgl. Abb. 9a, Ortskurven fiir Fall (1), (2) und (4)). 

Die Schnelle £, eilt der erregenden Kraft Fu 
bei tiefen Frequenzen um 90° voraus, im Gebiet 
der Resonanzstellen andert sich die Phase sehr 
stark zwischen 0 und —90° und bei hohen Fre- 
quenzen strebt die Schnelle einer Nacheilung bis 
Za OP ZzAle 


® 


© 


[/Fy| 
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SS 
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a 
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Abb. 9. (a) Ortskurve der relativen Schnelle an der Helm- 
holtz-Resonatoréffnung g L/Fu. 
(b) Ubertragungskurve |£,/F | =f(w/@m). 
(c) Ortskurve der relativen Bewegungsamplitude 
an der Helmholtz-Resonatoréffnung &,/F yy. 
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Der Absolutbetrag der Ortsvektoren dieser 
Kurve, multipliziert mit w/w, ergibt die gesuchte 
Kurve |é:/Fu| =f(w/wm) — vgl. Abb. 9b. Auch 
hier liefert der Fall (2) eine relativ ausgeglichene 
Ubertragung in einem groBen Frequenzbereich. 
Der Ubertragungsbereich ist um etwa 30% ge- 
geniiber einem einwelligen Schwingungsgebilde 
mit den gleichen Bestimmungselementen mm, 
sm und dy = 1 (vgl. Kurve 5, Abb. 9b) erweitert. 

Die mit —j@y/w multiplizierten Vektoren £,/Fu er- 
geben die Ortskurve der relativen Bewegung ¢,/Fyy an der 
Offnung des Helmholtz-Resonators — vgl. Abb. 9c. 

Bemerkenswert ist der Vergleich der Ampli- 
tudenabnahme oberhalb der Grenzfrequenz. 


Fall a 
Aus dem elektrischen Ersatzbild der Abb. 3 
]aBt sich ablesen, daB die Schnelle €y fiir hohe 


Frequenzen im wesentlichen durch den Schein- 
widerstand jwmm bestimmt ist: 


Ew 1 
as : 18 
(re). jomm v) 
Damit wird der relative Ausschlag 
Em Sail. 
= 1 
ae o2mm’ ah?) 
d.h. er fallt mit w? ab. 
Fall b 
Analog ist hier 
Eb jSL 
= = 20 
(Fe o-oo wmy(w*?mm — sz) 26) 
und da hier w?mm > Sz, ist: 
Eb SL 
= | 
laaley. w*mmymt’ ( ) 


d. h. die Amplitude fallt mit @‘. 

Im oberen Frequenzgebiet verhalt sich also der 
Ausschlag &y des Schwingungssystems wie der 
eines einwelligen Schwingungsgebildes, dagegen 
ist im anderen Fall der Helmholtz-Resonator als 
TiefpaBfilter in den Schallweg eingebaut und 
driickt damit die Amplitude &, starker herab. 

Kann man die Steife des Volumens Vr im Telephon 


nicht mehr vernachlassigen, d. h. wird sp < sy im Gegen- 
satz zu Gl.(8), so ist: 


Fall a sr 


W,=rpt+j (om ral (22) 
@) 

Unter Vernachliassigung der Realteile rj)y und ry, ergeben 

sich fiir den Betrag der Schnelle &), folgende Extremwerte: 


Eine Nullstelle fiir 


wo my—(sr+sr)=9, (23) 
zwei Unendlichkeitsstellen fiir 
(w*my—sm)[w*myp—(s_+ sr)—s_(w?mp—s7)|=9, (24) 


d. h. man erhalt hier zwei unwesentlich veranderte Reso- 
nanzstellen (vgl. dazu Gl.(13) und (15)). 
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Weiterhin ist der Ausschlag &y, bei tiefen Frequenzen 
wenig verkleinert, da hier 


(eal fej l SS 
EM /aiz> Oy igyes LSE 


7 Spu+ST 
ist. Bei hohen Frequenzen gilt hier G1.(19), so daB die 
Gesamtkurve der relativen Membranbewegung praktisch 
die gleiche wie die bei vernachlassigtem Hoérervolumen ist. 
Fall b 


Man hat zu setzen 


(25) 


+ 
Wu=rtj(omm— “MSP. (26) 
Fiir den Betrag der Schnelle Ey ergeben sich unter Vernach- 
lassigung der Realteile auch hier nur zwei Unendlichkeits- 


stellen, naémlich 
(27) 
d. h. die beiden Resonanzstellen sind wieder gegeniiber dem 
Fall vernachlassigbaren Hérervolumens unwesentlich ver- 
schoben (vgl. Gl.(15)). 

Bei tiefen Frequenzen ist der Ausschlag £; etwas mehr 
verkleinert als im Fall a, denn 


Ep dh 
Fee smtsp’ 
wihrend die Amplitude bei hohen Frequenzen wie im Fall 
des vernachliassigbaren Hérervolumens gemaf Gl.(21) mit 
w* fallt. Die Korrektur ist also in der Praxis so gering, da} 
man stets mit dem vereinfachten elektrischen Ersatzbild 


(Abb. 3) rechnen kann. 


[w*my—(smt s7)] (@2?m_—sy_)— w? my sp = 0, 


(28) 


3. Anwendung bei Mikrophonen 
Die fiir die Abstrahlung behandelten Fille 


lassen sich auch auf Aufnahmesysteme gleicher 
Bauart anwenden. Ein Kristallmikrophon soll als 
Beispiel dienen: Dieses enthalt einen an drei 
Ecken unterstiitzten Sattelbieger, der an der 
vierten Ecke mit einer Membran fest verbunden 
ist (Abb. 10). Bei einer solchen Anordnung liegen 
die beiden Federungen von Membran und Kristall- 
platte parallel (Abb. 11), so daf die wirksame 
Steife sy des Schwingungsgebildes 


(29) 
ist, worin sx die Kristallsteife sein soll. Als wirk- 
same Schwingmasse mp tritt die Summe der 


Sw = Sm-+ Sk 


aD 4 Sy iP 


ens) 


IO KIPK 
Abb. 10. Schematischer Aufbau eines Kristallmikrophons, 
dessen Membran mit einem Helmholtz-Resonator 
gekoppelt ist. 


Abb. 11. Mechanisches Ersatzbild des Mikrophons Abb. 10. 
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effektiv wirksamen Masse des Kristalls mx und 
die der Membran auf, also 


my = my + mx. (30) 


Und schlieBlich addieren sich die Reibungswider- 
stande der Membran ry und des Kristalls rx, so 
daB der wirksame Widerstand 
(31) 
wird. Mit diesen GréBen my, sy und rp ist aber 
das Gebilde grundsatzlich auf das eines unge- 
koppelten Schwingungsgebildes (z. B. wie an- 
fangs eine Membran) zuriickgefithrt. Auch hier 
ist wieder das Volumen des Gehduseinneren we- 
gen seiner sehr kleinen Steife vernachlassigt. 
Koppelt man an dieses Schwingungsgebilde 
einen Helmholtz-Resonator an und beschrankt 
sich bei der Berechnung auf einen Mikrophontyp, 
bei dem der Helmholtz-Resonator mit mr, sz und 
rr im Schallweg liegt, so greift die erregende 
Kraft Fy, an der Offnung des Helmholtz-Resona- 
tors an. Im mechanischen Bild (Abb. 11) liegt 
also diese Kraft zwischen der Masse my und der 
unendlich groB gedachten Gehausemasse. Die an 
dem Kristallsattelbieger auftretende Bewegungs- 
amplitude éx soll bestimmt werden. 


'w=—TkK+Irm 


thy Usy Cy ry m, n 


H 


muy mx 


Abb. 12. Elektrisches Ersatzbild des Mikrophons Abb. 10. 


In Abb. 12 ist das elektrische Ersatzbild 
2. Art aufgezeichnet. Hier liegt die Kraft Fy als 
elektrische Spannungsquelle in Serie zu den 
GréBen rz und my, wahrend die am Kristall auf- 
tretende Schnelle x dem Strom durch die 
Reihenschaltung der Tragerelemente der Membran 
mm, 1/sy und ry und der des Kristalls mx, 1/sx 
und rx entspricht. Gegeniiber der urspriinglichen 
Anordnung von Abb. 2 bzw. Abb. 3 ist also Ex 
und F,, vertauscht. Da man bekanntlich in sol- 
chen elektrischen Netzen Ursache und Wirkung 
vertauschen kann, ohne das Ergebnis zu andern, 
erhalt man also fiir diesen Fall die gleichen 
Frequenzkurven wie in b). 

(Eingegangen am 1. Februar 1953.) 
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VIBRATION OF PLATES COVERED WITH A DAMPING LAYER 


by A. VAN ITTERBEEK and H. MYNCKE 


Instituut voor Lage Temperaturen en Technische Physica, Leuven, Belgium 


Summary 


A study is made of the damping of steel plates covered with a thin layer of bitumen emulsion 
containing schist powder or similar products. This damping is characterized by the constant t 
which is the time corresponding to a decrease of the amplitude of the vibration to one half. After 
the description of the experimental set-up, results obtained with plates of various thicknesses 
and covered with different materials are given. The damping as a function of temperature, be- 
tween —20°C and + 80°C, is also determined. The check of the efficiency of the damping paint 
on a cement mill is also discussed. 


Sommaire 


On étudie l’amortissement de téles en acier enduites d’une mince couche de bitume, contenant 
du schiste moulu ou un produit similaire. Cet amortissement est caractérisé par la constante 7, 
c’est A dire le temps au bout duquel l’amplitude de la vibration est réduite de moitié. Aprés la 
description de l’appareillage, on donne les résultats obtenus avec des tdles de différentes 
épaisseurs et couvertes de différents matériaux. On détermine également l’amortissement en fonc- 
tion de la température entre —20° et 80°C. On a appliqué le procédé 4 un mélangeur de ciment 
enduit de l’un des produits en question: on a mesuré l’intensité et relevé le spectre sonore; les 
résultats semblent satisfaisants. 


Zusammenfassung 


Es wurde die Dimpfung von Stahlplatten untersucht, die mit einer diinnen Schicht Bitumen- 
lésung unter Zusatz von Schieferpulver oder ahnlichen Substanzen iiberzogen waren. Die Damp- 
fung wird durch die Gréfe t charakterisiert, das ist die Zeit, nach der die Amplitude auf die 
Halfte des Ausgangswertes abgeklungen ist. Nach der Beschreibung der Apparatur werden MeB- 
ergebnisse wiedergegeben, die an Platten verschiedener Dicke und mit Uberziigen aus ver- 
schiedenem Material gewonnen wurden. 

Ferner wurde die Temperaturabhangigkeit der Dampfung zwischen —20°--: + 80°C bestimmt. 

Die praktische Brauchbarkeit des Verfahrens wurde an einer Zementmischmaschine erprobt, 
die mit dem dampfenden Anstrich versehen wurde. Die Schallintensitét und das Schallspektrum 


wurden gemessen. Das Resultat schien zufriedenstellend zu sein. 


1. Introduction 


Because of its interest to the industry of motor- 
cars and shipbuilding we commenced two years 
ago in our laboratory a systematic study on the 
damping of steel plates, covered with a thin layer 
of a damping paint. This paint consists in general 
of a powder of very small particles of ground 


schist in a viscous emulsion of bitumen in water. - 


This mixture is spread over the plates by means 
of a spray-gun, the solvent evaporates after a few 
hours and on the plates there remains a thin 
layer of bitumen containing the fine powder. It is 
well known that even a very thin layer may pro- 
duce a very pronounced damping of the plate. We 
proposed to determine the influence of the thick- 
ness of the layer and also that of the temperature 
between + 80°C and —15°C. Recently there ap- 
peared a paper by OBERsT and FRANKENFELD [1], 
which treats a similar problem. However we are 
more interested in relatively thin layers, while 
the work of the previous authors is concerned 


more with layers deep compared to the thickness 
of the plates. 

The damping products, denominated by the 
symbols S,, 5,, S;, S,, and the steel plates were 
put at our disposal by the various industries. We 
have no indication how these paints were pre- 
pared, but an examination under the microscope 
revealed a particle size between 1 and 20 microns. 

In order to check the practical use of these 
paints, a cement mill was covered with product S, 
and measurements were made of the sound inten- 
sity and spectrum of the transmitted noise. These 
results are discussed under section 4. 


2. The experimental method 


The measurements were done on square plates 
of the following dimensions: 25 x 25 cm? or 
30 x 30 em?. The thickness varied between 1.5 
and 26 mm. To measure the damping we used 
the following mounting (see Fig. 1 and the photo- 
graphs Fig. 2, 3 and 4). 
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Fig. 1. Experimental set up used for the measurement of 
damping of steel plates. 


A 
Madd 


Fig. 2. The electronic equipment. 


The plate B, perforated in the centre, is fixed 
vertically by means of a brass rod (2 15 mm) 
which fits exactly in it. The rod rests with its 
extremities on two supports. 

The plate is excited by means of a wooden 
hammer C with a flexible haft, which is driven 
electromagnetically (key M). The vibrations of 
the plate are detected by means of the small 
coil E of an electromagnetic microphone. This 
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Fig. 4. Mounting of the plate in the refrigerated chamber. 


coil is placed at a distance of a few tenths of a 
millimeter in front of the edge of the plate, in an 
antinode of the vibration, in order to obtain the 
largest amplitude. This signal is filtered by the 
“Sound Analyzer” I (General Radio Co.) and 
after amplification connected to one pair of plates 
of a Cossor double beam oscillograph J. The 
“Sound Analyzer” is used in order to measure the 
damping of the fundamental vibration only. A 
study of the influence of a damping material on 
the harmonies of the fundamental frequency was 
also made. However the results have shown that 
the effect is rather small. These frequencies are 


ee 
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not so important, as in general the amplitude is 
smaller and the damping much greater. 

The measuring technique is as follows: 

In closing the contact M the hammer hits the 
plate. The signal picked up by the coil E appears 
then as a vertical strip of decreasing amplitude 
at the right hand side on the screen of the oscillo- 
graph. As the hammer is hitting the plate, con- 
tact D is closed and the electronic device K en- 
sures a linear horizontal displacement of the 
image to the left. In this way the damping curve 
will appear on the screen and can be photo- 
graphed. The single-stroke sweep given by K is 
obtained by the discharge of a condenser over a 
pentode tube RV2P800. The velocity of this 
single-stroke time base can be adapted to the 
damping by changing the condensers (switch Sc) 
and the variable resistor R. 

The time of the displacement of the spot from 
right to left over the screen can thus be varied 
between 12 seconds and a few hundredths of a 
second. 

The second beam of the oscillograph is con- 
nected to a timing device obtained by the dis- 
charge of a condenser over a neon tube. The fre- 
quency of this time indicator is changed by 
means of switch Sr. The calibration is done by 
means of a chronometer for the three first posi- 
tions and photographically with reference to the 
power supply frequency for the following ones. 
The five positions normally used give the follow- 
ing values for the relaxation period of the neon 
tube system: 


Position number: 1 2 3 4 5 
Relaxation ; ‘ 
period: 1.84s 1.18s 0.51s 0.17s 0.045 


After arrangement of sweep velocity and time 
indicator frequency the damping curve is photo- 
graphed with a Cossor Camera. 


Fig. 5. Picture obtained on the Cossor screen. 
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Then the spot is brought to the right by means 
of key L so that the next picture can be taken. 

From this photograph the damping constant t 
is then determined. This means the time Tt cor- 
responding to a decrease of the amplitude A to 
one half, if it is supposed that the decrease is 
logarithmic. t is calculated from the slope of 
log A as a function of time by means of the 
equation: 

0.301 
A (log A) /At’ 


To measure the influence of temperature on the 
damping between room temperature and 80°C 
and also between room temperature and —15°C, 
the plates with the excitation system are mounted 
within two thermostats. The first is kept constant 
at temperatures between 20°C and 80°C, the 
second is a large scale household refrigerator. 

To measure the temperature a small chromel- 
alumel thermocouple F is fixed to the plate and 
connected to a galvanometer H as indicated in 
Fig. 1. 

All the controls for taking the pictures are 
operated from outside the experimental chamber. 


3. Experimental results 
a) Comparison between different products 


In Fig. 6a-k are graphically drawn the ex- 
perimental results obtained for different kinds of 
products: S,, S,, S, and S,. 

The products are spread out on a plate of 
1.5mm thick (250x250 mm?) and on steel 
plates used in shipbuilding, put at our disposal 
by John Cockerill Co., of 300 x 300 mm? and the 
following thickness: 3; 4; 5.3; 6.3; 7.8; 8.2; 9.1; 
10.5; 20.0 and 26.0 mm. We measured the de- 
crease of the constant t as a function of the 
weight of the material spread out on the total 
surface of the plate. In the graphs we give the 
number of grams of the wet material. The effect 
of drying decreases approximately this weight 
with the following percentages: 


material: oH Se 5. D4 
20% AL On S4 er, | 34% 


decrease in weight: 


The fundamental frequency of the plate is indi- 
cated on each figure. The point on the ordinate 
axis corresponds to the bare plate. Taking into 
account the average density of the materials 
used, it can be said that the thickness of the layer, 
spread over an area of 900 cm?, is approximately 
0.7 mm for 100 g of the wet material. 

From the various experimental results, the 
following conclusions can be drawn: 
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250% 2501.5 mm3: 88 c/s: A= S1; O = So 
e= $3; A= Sa 


300 *300x3mm3- 68 c/s: O = Sp: @ = S3 


300x300x4mm3; 143cls ; O = So; @= Sy 


0 40 80 120 160 


300x300x 6.3mm3; 235¢c/s; © = Sp 


300x300x7.8mm3; 270 c/s, © = So 


P.LIENARD: LA RESONANCE D’UN OSCILLATEUR 


ACUSTICA 
Vol. 3 (1953) 


300 *300x9.1 mm3; 340c/s; O = So 


[ 


300x 300x10.5 mm3, 379 c/s ,O =S9 


rina 


0 20 40 60 80 100 g 
5s 2] 

[ 300%300x 20 mm3; 690 c/s; & = Su 
4 


| 


0 
3s 
| | | | 300%300%26mm3; 950 c/s; A= Sy 
3 =rHl T 7 jase | 
2 | | Le} 1 
©, ® 
0 20 40 60 80 1004 0 40 80 120 160 2004 
—= mass 


Fig. 6. Damping for different masses of material spread on the plates. 


Large differences appear between the de- 
creases obtained for the different products. This 
difference is more pronounced as the thickness of 
the plate increases. For the good products a 
sharp reduction of t is obtained for a covering 
weight smaller than 25 g. For the plates of which 
the thickness is between 6 and 1] mm, the maxi- 
mum reduction of t obtained is about 50%, 
while for the thicker plates this reduction is much 
larger, but the fundamental frequency is also 


higher. This agrees with our practical measure- 
ments on the damping of a drum of a cement mill 
as reported under 4. 


b) Measurements as a function of temperature 


The experimental results obtained for the dif- 
ferent plates covered or not covered are indicated 
respectively in Fig. 7a-c. For the bare plate we 
observe a strong increase of t as the temperature 
decreases. 
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300*300x3mm3; 4 = bare 
© =53(107g) 


300x300x4mm3; A = bare 
©=Sy (51g); @= $3(83g) 


300x300x8.2mm3; a = bare 
©=$7(73g); © = $3(1I1g) 


-20 0 20 40 60 
—~ temperature 


Fig. 7. 


Variation of damping with temperature 


80 °C 


Further we can observe, that from about 40°C 
the effect of covering the plate is practically nil; 
this is probably due to the fact that the paint 
becomes soft in the neighbourhood of this tem- 
perature. However, between room temperature 
and —20°C the effect is very large, so that we can 
conclude that this damping paint will be more 
efficient in the winter than in the summer. 


4. Practical check of the efficiency of the damping 
paint on a cement mill 


As is well known, the sound level produced by 
a cement mill is very high (about 100 dB). In 
order to check the practical use of these damping 
paints we asked a firm in the cement industry to 
cover one of their mills with a layer of one of the 
paints, which we estimated from our previous 
measurements to produce the largest effect, 
namely S,. This cement mill has a length of 6 m 
and a diameter of 3.8 m. The thickness of the 
steel plate is 26 mm. 

We measured, by means of a sound level meter 
of the General Radio Co., the sound intensity in 
the neighbourhood of the uncovered mill. The 
analysis of the sound by means of a Sound Ana- 
lyzer of the General Radio Co. was also carried 
out. The drum of the mill was then covered with 
a layer of 4 kg/m? of wet material, which cor- 
responds to a thickness of the dry layer of about 
2.9 mm. 


100% 
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Fig. 8. Sound spectrum produced by a cement mill. 


The sound produced was analysed again. The 
spectra obtained are represented in Fig. 8. The 
relative amplitudes are expressed in percentages 
against frequency, where 100% is taken for the 
loudest component. The full line corresponds to 
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the uncovered drum, while the dashed curve 
refers to the covered mill. For the sound level, 
we found 95 dB in the first case, while in the 
second one a decrease of 5.7 dB was obtained. 
Moreover we observe that the change in the 
spectrum is rather important. Also by simple 
audition, this effect could easily be observed. 
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REPRESENTATION GRAPHIQUE DE LA RESONANCE 
D’UN OSCILLATEUR, UTILISABLE POUR L7ETUDE DE 
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Sommaire 


Si l’on prend pour paramétres l’accélération et l’amplitude en coordonnées homogénes d’un 
oscillateur linéaire 4 un seul degré de liberté, la courbe de résonance est une ellipse qui permet de 
retrouver les propriétés fondamentales classiques de la résonance. 

Si l’oscillateur n’est pas linéaire, la figure de résonance est un lobe plus ou moins déformé de 
Vellipse. On met ainsi en évidence l’influence de certains paramétres de non linéarité, méme trés 
petits, et on indique des applications expérimentales possibles. 


Zusammenfassung 


Nimmt man als Parameter die Beschleunigung und die Amplitude (in homogenen Koordinaten) 
eines linearen Oszillators mit einem Freiheitsgrad, so ist die Resonanzkurve eine Ellipse, mit 
Hilfe derer man die klassischen Resonanzeigenschaften auffinden kann. 

Ist der Oszillator nicht linear, so wird die Resonanzellipse mehr oder weniger verzerrt zu einer 
allgemeineren Schleifenform. Man erkennt daraus leicht den Einfluss gewisser Parameter der 
Nichtlinearitat, auch wenn diese sehr klein sind; mégliche experimentelle Anwendungen werden 


angegeben. 


Summary 


If acceleration and amplitude are taken as the parameters (in homogeneous co-ordinates) of 
a linear oscillator with only one degree of freedom, the resonance curve is an ellipse by means of 
which the fundamental properties of resonance can be derived. 

If the oscillator is not linear, the resonance ellipse is distorted and lobes appear. Thus, the 
influence of certain—even very small—parameters of non-linearity is pointed out; possible 


experimental applications are mentioned. 


1. Introduction 


L’oscillateur classique, constitué d’une masse 
m, d’une raideur K et d’un frottement propor- 
tionnel a la vitesse, de coefficient r, excité en 
vibrations forcées par une force sinusoidale de 
valeur maximum P, est parfaitement connu 
théoriquement car il conduit a l’équation facile- 
ment intégrable: 


* Réimprimé avec quelques changements non essentiels 


de «La Recherche Aéronautique» 29 [1952], 45. 


mX +rX + KX = Psinat 


ou, avec variables et constantes homogénes, en 
posant: 


2 


K/m= M0, 


P= KX,=mX,o 


r/2m = &Wo, 


et en remplagant X/X, par x: 
7 5 2 ee 
x + 2ew yx + Wox = wo sin wt 


(1.1) 


dont la solution est: 
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x =asin(wt — q), 
2\ 2 27 —% 
es I(1- <3) +4ee 
Mo ay) (1.2) 
a) gives 
teg=—2e i : 
“A al na) 


Les courbes représentant l’amplitude relative a 
en fonction de l’écart des fréquences w/w, sont 
les courbes de résonance bien connues, dont la 
forme varie suivant les valeurs de l’amortisse- 
ment €. 

On sait que l’amplitude a passe par un maxi- 
mum pour une certaine valeur de la fréquence 
excitatrice, inférieure a la fréquence propre du 
systéme non amorti, et que le maximum de 
vitesse et le maximum d’accélération 
s’obtiennent pour des fréquences différentes 
(Fig. 2). 

Pratiquement, ces diverses fréquences sont trés 
proches les unes des autres tant que l’amortisse- 
ment n’est pas trés élevé, et on confond ces trois 
maxima, ainsi que l’égalité de la fréquence excita- 
trice et de la fréquence propre, sous le nom de 
résonance. 

La représentation ci-dessous nous permettra, 
entre autres, de mieux distinguer ces divers phé- 
nomeénes. Le rapport a w, de Pécart Aw des deux 
fréquences pour lesquelles l’amplitude est égale 
a la moitié de l’amplitude maximum est tel que 


do V3 AW 
= ~a 


iN 
Wy 2x 3) 
AW/W étant la perte relative d’énergie 
par cycle de lVoscillateur, ou capacité 
d’amortissement, égale par ailleurs au double 
du décrément logarithmique 0: 


AW 
WwW = 20 = 47 é. (1.4) 
Le rapport 
Le fete 
— eae (1.5) 


est appelé amplification a la résonance de l’oscil- 
lateur, qui est aussi égal au coefficient de sur- 
tension. 

Ces diverses quantités peuvent servir a carac- 
tériser l’amortissement [3]. Mais on oublie sou- 
vent que si certains oscillateurs (circuits élec- 
triques, masses suspendues élastiquement et 
amorties dans un liquide, etc.) réalisent avec une 
approximation suffisante les conditions de validité 
de ’équation (1), la plupart des oscillateurs 
mécaniques dont l’amortissement est da 
a un frottement intérieur ne satisfont 
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pas a l’équation linéaire, parce que le frot- 
tement n’est généralement pas proportionnel a la 
vitesse. C’est le cas, par exemple, de parois 
minces, plaques, diaphragmes ou cloisons qui 
vibrent sous l’influence d’une onde sonore, avec 
une amplitude qui dépend des fréquences propres 
de ces parois, et aussi de leur amortissement. 

Dans ces oscillateurs, le décrément sera fonc- 
tion de l’amplitude et non plus une constante. 

Enfin, la force de rappel élastique peut étre 
aussi un terme non linéaire, ce qui suppose un 
écart a la loi de HooKE. 

Les équations des systémes non linéaires sont 
assez peu connues, étant souvent non-intégrables, 
sauf cas particuliers. Par suite on cherchera le 
plus souvent 4a faire rentrer «de force» le fonc- 
tionnement d’un oscillateur pratique dans le 
cadre des équations linéaires, avec dissipation 
proportionnelle a la vitesse et force de rappel pro- 
portionnelle a l’amplitude, ce qui équivaut a faire 
une approximation parfois dangereuse. 

C’est ainsi que dans |’étude de l’amortissement 
d’une barre d’acier, dont le décrément est cons- 
tant avec une trés bonne approximation, on 
reléve cependant un rapport de l’ordre de 1 a 3 
ou 4 entre le décrément mesuré directement et 
celui calculé a partir de la courbe de résonance. 
Cette observation nous a conduit a étudier de 
plus prés cette résonance, et a mettre en évidence 
un terme de non linéarité du rappel élastique de 
Vacier, dont Vinfluence est suffisante pour ex- 
pliquer cette divergence. 

L’étude théorique exposée ci-dessous a pour but 
de donner une représentation pratique du phéno- 
méne de résonance, dans le cas de systémes non 
linéaires, quand on ne cherche pas a connaitre la 
forme méme de la vibration, mais seulement les 
valeurs maxima de |’amplitude ou des efforts. 

Elle a pour base l’observation faite ci-dessus 
et une méthode d’approximation indiquée par 
TIMOSHENKO [1]. 


2. Représentation de la résonance dans le cas 
de systéme linéaire 


2.1. Transformation fondamentale 


La représentation graphique exposée ci-dessous 
permet de mieux figurer les phénoménes de réso- 
nance méme dans le cas bien connu d’un systéme 
linéaire. L’application aux autres systémes n’en 
sera qu’une extrapolation, avec une approxima- 
tion a la base évitant l’intégration. 

Nous écrirons d’abord l’équation (1.1) sous la 
forme: 


“ bd ran 2 jos 
x+ ew yx + wox = aoe” (2.1) 
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On sait que lintégration de cette équation 
revient a écrire: 
P= IDX, KL — OX 
Cane a ; 
donc: —(—)'xt2+2je—az=e™. (2.2) 
Wo Wo 


Explicitant le module de x par x= |x|- eliot—o), 


divisons d’abord (2.2) par e’”. Egalons alors 
les parties réelles et imaginaires des deux 
membres, et élevons au carré chaque équation. 
En faisant la somme, il vient une relation simple 
entre | x| et (w/ap)”. 

Nous lécrirons en convenant que x désigne 
maintenant la valeur absolue de l’amplitude, 
pour simplifier l’écriture, et en posant de plus 


(o/a)?=y|%5 
(x—y)?+ 4e2xy —1=—0. 


on obtient (2.3) 
Cette derniére transformation revient a écrire 
l’équation conjuguée de (2.2) et a faire le produit 
des deux équations. 

La méme équation (2.3) s’obtient d’ailleurs 
par léquation (1.2) en remplacant a par x et 
(w/o)? par y/x. 

Done dans (2.3), x représente la valeur ma- 
ximum de l’amplitude, et le rapport y/x=tg0 
est le carré du rapport de la fréquence d’excita- 
tion a la fréquence propre du systéme non amorti. 
D’autre part, on peut écrire: 
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y= xw?/oo 


c’est donc le rapport de l’accélération xm? a la 
limite w; de cette valeur pour w augmentant 
indéfiniment (voir éq. 1.2) ce qu’on pourrait ap- 
peler la surtension instantanée ou l’accé- 
lération relative (a celle imprimée par une 
pulsation infinie), x étant l’amplification en am- 
plitude, ou amplitude relative (a celle impri- 
mée par une pulsation nulle). 

L’équation (2.3) définit un faisceau de coniques 
centrées sur l’origine (Fig. 1). Les points de base 
ont pour coordonnées: x= +1, y=+1. 

Les axes de symétrie sont: y= +x. 

Pour «=0 (aucun amortissement) la conique 
se décompose en deux droites y=x+ 1. 

Pour ¢ compris entre 0 et 1/J/2, c’est une ellipse 
dont le grand axe est y= x. Sie=1/J/2, amortisse- 
ment pour lequel l’amplification a la résonance 
est égale a 1, l’ellipse devient le cercle x*+ y?=1. 

Pour ¢ compris entre 1/2 et 1, l’ellipse a son 


grand axe suivant y=—x. 
Pour ¢=1 (amortissement critique), elle se 
décompose en deux droites y=—x +1. 


Pour ¢> 1, la conique est une hyperbole. 

Le lieu des amplitudes maxima (dy/dx= 0) est 
Vhyperbole équilatére x?—y?—1=0, et le lieu 
des surtensions maxima dy/dx=0 est l’autre 
hyperbole y?—«?—1=—0. La Figure 1 représente 


le faisceau de coniques pour diverses valeurs de e. 


Fig. 1. Représentation des ellipses de résonance d’un systéme linéaire pour diverses valeurs de l’amortissement; 


x 


Si 
o/@» 


Hoi il 


amplitude relative, 
accélération relative, 
pulsation relative. 


Graduation par méthode graphique de l’axe de pulsations. 
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2.2. Etude de la résonance au moyen de cette repré- 
sentation 


Il est facile maintenant d’interpréter cette 
figure et de passer graphiquement 4 la représen- 
tation classique des courbes de résonance (Fig. 2). 


Fig. 2. Construction d’une courbe de résonance classique 
x = f(w/m) a partir de Vellipse de résonance 
(exemple ¢ = 0,2). 


L’amplitude est donnée directement en va- 
leur relative par l’abscisse x. Une demi-droite 
issue de l’origine, a pour pente tg = w?/a. 
Done seul le premier quadrant est a considérer 
physiquement, et la pente de cette droite est 
fonction de la fréquence. L’axe y= x correspond 
a Végalité de la fréquence d’excitation et de la 
fréquence propre. Dans le cas ot «= 0, c’est-a-dire 
quand l’oscillateur n’est pas amorti, l’amplitude 
croit depuis x= 1 pour la fréquence nulle, jusqu’a 
Vinfini pour w=q@, et décroit ensuite jusqu’a 
zéro. On peut graduer en valeur de w un cercle 
centré en O, et construire par points la courbe de 
résonance. 

Pour 0<e<1///2, la courbe de résonance 
passe par un maximum non infini. 

On distingue nettement les trois points: maxi- 
mum d’amplitude, égalité des fréquences, maxi- 
mum d’accélération. C’est Pégalité des fréquences 
qui est couramment appelée résonance (résonance 
de vitesse); l’amplitude maximum en est nette- 
ment différente ainsi que l’accélération maximum 
(résonance d’accélération). Le maximum d’am- 
plitude est obtenu par la pulsation «, telle que 
Véquation (2.3) ait une racine double en y, soit: 


y= Wy) (1—2e?)”. (2.4) 


De méme le maximum d’accélération est 


obtenu pour 


Wg= Wy (1—2e?)—”. (2.5) 
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Ces résultats sont classiques, mais retrouvés ici 
trés rapidement. Si l’on coupe l’ellipse par une 
droite x= const., l’examen du diamétre conjugué 
de cette direction conduit 4 la relation suivante 
entre les pulsations w,, ; correspondant A la 
méme amplitude x (Fig. 3): 


wo, + ow =203= 205 (1 —2¢*), 


(2.6) 


ce qui peut s’exprimer ainsi: 

La somme des carrés des pulsations 
donnant la méme amplitude est cons- 
tante et égale A 2m, (1—2e?) quelle que 
soit l’amplitude. 


@; 
y ae 
w, 
W@W) 
M 
1 A 


——— 1X 


Fig. 3. Etablissement graphique de la relation (2.6) 


MI = MA + ATI, IN = TA + AN, MI =T 


En considérant le diamétre conjugué de la 
direction y=const., on obtient de méme: 


Oo +o," =20,=20,(1—2e*)—.. (2.7) 


La somme des carrés des pulsations 
donnant la méme accélération est cons- 
tante et égale a 2w,(1—2e?)—. 


2.3. Etude de la capacité d’amortissement 


Plus particuliérement, si on désigne par w, et 
w; les pulsations correspondant A |’amplitude 
égale a la moitié de l’amplitude maximum (pour 
laquelle la pulsation est w,) le calcul de la diffé- 
rence @; —@, = Aw conduit, en prenant pour 
(1—e?)” un développement limité, 4 l’expression: 


Aw |= 2¢ V3 + 422 V3 + ke. (2.8) 


Expérimentalement, il est impossible d’utiliser 
la pulsation w, du systéme non amorti. On em- 
ploiera, soit la pseudo-pulsation wy, soitla pulsa- 
tion w, correspondant au maximum d’amplitude. 
Comme @)=)(1—e)% et w= @)(1—2e?)% 
on obtient: 


Aw|a, = 2 y3(1 + 
Aww, = 2¢ V3 (1+ 322+ ke*). 


e?+ ke®), 


(2.9) 
(2.10) 
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Ces trois expressions (2.8, 2.9, 2.10) se réduisent 
a la formule unique 


Aw|w, =2eV3 
quand « est petit et que l’on confond les diverses 
pulsations de résonance. 

La relation indiquée plus haut (1.4) permet de 
relier la capacité d’amortissement AW/W au 
rapport Aw/m,) mesuré expérimentalement, par 
cette équation (2.11). 

L’expression (1.4) n’est, elle aussi, qu’approxi- 
mative. Son calcul se fait, a partir de l’équation 
(1.1), sans deuxiéme membre, en évaluant l’énergie 
de Voscillateur égale 4 la somme de l’énergie po- 
tentielle et de l’énergie cinétique, dont on prend 
la valeur moyenne, puis la perte d’énergie, égale 
au produit par la vitesse de la force de frottement 


(2.11) 


rx, dont on prend également la valeur moyenne 
par cycle. 
On obtient ainsi pour la capacité d’amortis- 


sement (2.12) 


—— = Agre(1—e?) * =4ate Tee ee 
W 2 8 


en remplacement de |’équation 1.4. 

Le calcul plus précis de AW/W se fait a partir 
des équations (2.8)... (2.12), par approximations 
successives. 

On trouve, en définitive, pour expression de la 
capacité d’amortissement, suivant que l’on 
rapporte la largeur Aw de la courbe de résonance 
a la pulsation m,) du systéme non amorti, a la 
pseudo-pulsation w, ou a la pulsation du maxi- 
mum d’amplitude @,: 


Aiea 2F yelp I 0,095 (S) +. (2.13) 
W V3 ®o o 


2 
Aue BAeNeD k 0,136 ea tel, (2.14) 


Wie aS Sai Wo 


S 


iJ 


AU. Bw 0,178 (a) +4 (2.15) 
W V3 Wo 


Ces trois expressions se raménent usuellement a 


Ay apes 
W ." V3 Wo 
telle que nous l’avons employée pour la détermi- 


nation expérimentale de la capacité d’amortisse- 
ment, tant que celui-ci n’est pas trés grand. 


(1.3) 


2.4, Détermination des fréquences propres 


Pour les oscillateurs linéaires, la courbe de 
résonance relevée expérimentalement permet de 
construire l’ellipse de résonance (2.3). En effet, 
on peut déja caleuler ¢ par ’équation (2.10), ou 
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bien construire par point la courbe de l’accéléra- 
tion en fonction de la fréquence (pour chaque 
pulsation m, amplitude étant a, l’accélération 
est aw?) et utiliser les relations (2.6) et (2.7) pour 
deux couples de points de méme amplitude ou de 
méme accélération. 


On déduit en effet de (2.6) et (2.7): 


Dwr 5 
5 = (1—2e?)? et Loi La; = 4a. 


= 2.16 
So! (2.16) 


On obtient ainsi, encore la valeur de «¢, et celle 
de w,. Cette derniére permet de graduer en va- 
leurs de w/w, V'axe des abscisses de la courbe de 
résonance expérimentale, et de passer ensuite a 
la construction de l’ellipse par l’opération inverse 
de celle indiquée plus haut. On peut encore cons- 
truire graphiquement l’ellipse connaissant les 
deux points de base x=1, y=0 et x =O} 92 
et la tangente x= max. (il faut pour cela con- 
naitre la déformation statique X,). Il est bien 
évident que cette construction est pratiquement 
impossible si ¢ est trés petit, car ellipse est alors 
extrémement allongée. Dans ce cas, on ne fait 
d’ailleurs pas une grande erreur en prenant pour 
@, la pulsation du maximum d’amplitude. La 
construction indiquée présente un intérét dans 
le cas d’amortissements importants, et pour véri- 
fier si un oscillateur peut étre considéré comme 
linéaire. Elle permet également (soit a partir de 
l’équation, soit a partir d’une courbe expérimen- 
tale) de déterminer le maximum de vitesse ou le 
maximum d’accélération, parfois plus utiles a 
connaitre que le maximum d’amplitude. 


3. Systéme non linéaire avec amortissement 
visqueux 


3.1. Equation générale, principe de l’approximation 


Supposons d’abord a notre systéme une force 
de rappel proportionnelle a F(X), fonction quel- 
conque de l’amplitude X, l’amortissement restant 
proportionnel a la vitesse. 

L’équation différentielle est alors: 


mX +rX + KF(X)= Psinot. 
Nous poserons comme précédemment: 
K/m=a}, r/2m=ew), P= KX,=mX,on 


et nous rendrons les variables homogénes en 
posant: X/X,—-x. Il est a4 remarquer que nous 
rapportons maintenant les amplitudes a la défor- 
mation statique de l’oscillateur dont on aurait 
annulé les termes non linéaires (la déformation 
statique de Voscillateur réel correspondant a 


P=KF(X,). 
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D’autre part, il est bien évident que les coefli- 
cients apparaissant dans F(X) ne pourront pas 
étre tous sans dimension, mais que certains d’entre 
eux dépendront de l’amplitude X,, c’est-a-dire 
de la force appliquée. En revenant a des coeffi- 
cients m, z, K constants, l’ellipse de résonance 
précédente sera remplacée par une famille de 
courbes ayant pour paramétre la force excitatrice. 
Avec le changement de variable indiqué plus 
haut, et en posant de plus (1/X,) F(X)=f(x), il 


vient: 
x+2ewyx+o, f(x) =a, e™. (3.1) 


La solution exacte d’une pareille équation est 
inconnue. Des solutions approchées pour |’équa- 
tion sans second membre et e=0 (oscillations 
libres non amorties) sont indiquées, par méthode 
graphique ou par approximations successives, 
par TimosHEenKo [1]. L’oscillation n’est plus 
sinusoidale mais comporte des harmoniques. 

La méthode de résolution approchée consistera 
A essayer une solution de la forme x=ae™ et a 
déduire de (3.1) une expression de x en fonction 
de wt, d’ow on tire une expression meilleure de x 
par double intégration. 

Sila valeur maximum de I|’amplitude est 
seule recherchée, il n’est pas nécessaire de con- 
naitre les coefficients des divers harmoniques. On 
peut alors se contenter, comme l’a_ exposé 
TIMOSHENKO, dans un cas simple, de vérifier 
Véquation différentielle du systéme avec x= 
acoswt quand x prend sa valeur maximum, c’est- 
a-dire pour cosmt= 1. 

Avec les notations précédentes, ceci revient a 
satisfaire V’équation (3.1) en posant x—jox, 
x=—w?x, comme on I’a fait pour l’intégration 
exacte de l’équation linéaire (1.1). x désigne 
maintenant la valeur maximum del’ampli- 
tude, et on posera encore (w/a .)?=y/x, aprés 
avoir tout divisé par w5, c’est-a-dire qu’on rem- 
place x par —y, x par j(xy)” et w, par 1. 

En effectuant ces transformations sur ]’équa- 
tion (3.1) et en écrivant l’équation conjuguée il 
vient: 


f(x) —y — 2je(xy)* =e™, 
f(x) —y + 2je(xy)* =e" 
et par produit membre 4 membre: 
[f (x) —yP + 4e?xy = 1. 


Nous désignerons cette courbe, qui remplace 
ellipse précédente, par «figure de réso- 
nance», afin de la distinguer de la «courbe de 
résonance» usuelle tracée dans les axes (x, «). 


(3.2) 
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3.2. Utilisation de cette représentation (Fig. 4) 


L’équation (3.2) remplace celle de l’ellipse 
précédente (2.3). Les axes x, y et la pente 
y/x=tg 0 ont les mémes significations. S’il n’y a 
pas d’amortissement (¢= 0), la courbe se décom- 
pose en deux autres ayant en général deux 


branches infinies, y—f (x)= +1. 


Cy-flx) +4621 =1 


Fig. 4. Figure de résonance d’un systéme non linéaire 
(mémes axes que Fig. 1). 


Pour ¢> 0 la courbe (que nous limiterons en- 
core au quadrant x > 0, y > 0) est coupée en deux 
points par une paralléle 4 Oy (points de méme 
amplitude) et le lieu du milieu de ces deux points 
est la courbe y=f(x)—2e?x qui remplace le dia- 


300 


N 
S 


_--}------4------4+ 
ar 1 1 


Amplitude 


8 


fe es mete 


0 


1 3 5 
— Graduation arbitraire 
fonctionde wkx 


Fig. 5. Courbe de résonance expérimentale obtenue avec 
une barre d’acier trés peu amortie, tenue en ses 
neuds de vibration, correspondant a une figure de 
résonance analogue au cas de la Fig. 4. 
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métre conjugué de x=const. par rapport aux 
ellipses précédentes, obtenues pour f(x) =x, mais 
une demi-droite y/x—w?/w, peut couper la 
courbe en plusieurs points, ce qui correspond sur 
la courbe classique de résonance, a des sauts 
d’amplitude (v. Fig. 5 et [2] pour une courbe ex- 
périmentale de ce type). 


4, Etude d’équations particuliéres 
4.1. Systéme a@ rappel élastique de degré impair 


Nous étudierons plus 4 fond un type d’oscilla- 
teur dont se rapprochent plusieurs oscillateurs 
mécaniques réels, et présentant les caractéris- 
tiques suivantes: 


— la force de rappel élastique comporte un terme 
linéaire, et des termes de degré impair, dont le 
premier seul sera conservé, et considéré comme 
un terme correctif, c’est-a-dire de coefficient 
trés petit. Un tel systéme est symétrique par 
rapport a4 sa position d’équilibre (la force 
change de signe avec |’élongation). 

— lamortissement comporte un terme «vis- 
queux», proportionnel a la vitesse, et un 
terme proportionnel au produit de la vitesse 
par l’élongation. 

L’équation différentielle correspondante pour 

le systéme en vibration forcée sera: (4.1) 

t 


23 ae : - 5 
x + ZewMyx + Hox + 2yM_xx + 10) xX°> = we” 


par un changement de variables analogue a celui 
de (3.1), y et ~ étant respectivement proportion- 
nels aux puissances | et 2 de X,. La transforma- 
tion précédente conduit 4 l’équation: 


(x + xx? — y)? + day (e + yx)? = 1. 


6 ie 


y-0 Ye- he 


(4.2) 


Fig. 6. Figures de résonance — Equation (4.2) pour « = 0, 
6 = Olety=.07 0,012.0 reti0.5. 
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Les figures correspondantes et les courbes de 
résonance qui s’en déduisent graphiquement sont 
construites en Fig.6 avec les valeurs: 1=0; 
e= 0,1; y=0; 0,01; 0,1 et 0,5; et en Fig. 7 avee 
les. valeurs: e=0,1; y=0;.0=—0;,.0,0]l5,0 50a 
—0,1; —0,01. 


8 x 


Fig. 7. Figures de résonance — Equation (4.2) pour e= 0,1, 
y = O'et an = 0; 0,01 eb OF) 


I] est 4 remarquer que dans le cas d’amortisse- 
ment faible, une valeur méme peu élevée du 
coefficient « du terme non linéaire entraine un 
écart important entre la fréquence propre w, du 
systéme linéaire et la fréquence de résonance du 
systéme non linéaire. 

Dans le cas de y=0, on précisera au chapitre 5 
la forme de ces courbes. 


4.2. Oscillateur de relaxation de VAN DER PoL 
L’équation classique de VAN DER POL s’écrit: 
- : 
x — 2ewy (1 — Bx”) x + wx = 0 


que l’on peut écrire encore, en faisant figurer la 
force extérieure excitatrice: 


x—2ew)x + 2Bayextx+o,x=o,e”. (4.3) 
En effectuant la transformation précédente 
x>—y, x —> j(xy)”, W—>1 
on obtient: 


(x—y)?+ 4e2xy(1—Bx?)?=1. (4.4) 
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L’équation obtenue est indifférente au signe de 
é, c’est-a-dire qu’une attaque donnée produira la 
méme amplitude sur le systéme amorti ou sur le 
systéme auto-oscillant 4 amortissement négatif 
de méme valeur absolue. 

Les courbes obtenues sont analogues a celles 
de la Fig. 6 (cas précédent pour «= 0). 

On voit nettement sur l’équation (4.4) comment 
se stabilisera l’amplitude, du fait du coefficient /: 
pour x croissant, 1— fx? arrivera a la valeur zéro, 
annulant ainsi le terme d’amortissement. L’oscil- 
lateur vibre donc sur sa fréquence propre comme 
un oscillateur linéaire sans amortissement (ellipse 
de résonance décomposée, x —y = + 1), avec cette 
différence que l’amplitude n’est pas infinie mais 
déterminée par Bx?= 1. Pour f petit, cette ampli- 
tude d’oscillation stable sera d’ailleurs inférieure 
a 1. Le point de fonctionnement sur la figure de 
résonance se placera sur y=x+ 1. 


4.3. Résolution plus exacte de l’équation (4.1), pre- 
miére approximation 


La transformation précédente 


x—>—y, x—>j(xy)” 

revient a ne considérer que |’amplitude du sys- 
téme correspondant a la composante fondamen- 
tale X=ae!“ en phase avec la force excitatrice, 
en négligeant, d’une part, les harmoniques suc- 
cessifs, d’autre part, la composante en quadrature 
avec l’excitation. Pour une approximation meil- 
leure, et suivant la méthode indiquée par T1mo- 
SHENKO on peut chercher 4 satisfaire |’équation 
par une fonction de la forme: 


X=aXysinwt+ bX,coswt (4.5) 


en mettant cette équation (4.1) sous la forme 
x=f(x,x) en y remplacant x et x par les valeurs 
tirées de (4.5), ce qui conduit en intégrant terme 
a terme, 4 une nouvelle forme pour x. II serait 
possible de réitérer cette approximation, mais le 
calcul devient vite inextricable. Pour la premiére 
approximation seulement, il ne présente aucune 
difficulté particuliére, mais il est déja trés fasti- 
dieux. 

On obtient, en définitive, pour la nouvelle 
valeur de x, plus exacte, en regroupant les termes 
de méme pulsation, sous la forme d’un dévelop- 
pement en série de Fourier: 
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On voit que la racine de (4.1) ne peut comporter 
de composantes harmoniques paires si y= 0. Les 
coefficients des composantes successives sont 
d’autant plus petits qu’on s’éléve plus dans l’ordre 
des harmoniques. 

L’expérience montre que, pour des systémes 
non linéaires peu amortis, la courbe de résonance 
présente une distorsion trés nette, conduisant 
méme 4 des sauts d’amplitude, comme sur la 
Figure 5, sans que l’oscillation comporte un taux 
d’harmonique notable. C’est-a-dire que l’on peut 
considérer déja comme une assez bonne approxi- 
mation |’équation (4.6) réduite au premier har- 
monique, donc connaitre la valeur maximum de 
Vamplitude en égalant respectivement a a et b les 
coefficients de sinwt et coswt, les autres termes 
étant négligés. 

On calculera ainsi, d’abord a et 5b en fonction 
de a?+ 6b? puis cette amplitude maximum égale a 
(a2+ 62)”. Utilisant les variables précédentes: 


(a2 + B)* =x, wolo® = x/y 


on trouve facilement: 


As?xy = 1] (4.7) 
équation analogue a (4.2) (on a négligé les termes 
en y) dans laquelle le calcul plus précis montre 
qu’il faut remplacer « par 3x/4. 


4.4, Deuxiéme approximation 


La résolution suivant l’équation (4.6) fait ap- 
paraitre des valeurs approchées des composantes 
harmoniques 2 et 3, composantes que nous consi- 
dérerons comme des termes correctifs 4 introduire 
dans |’équation (4.7) d’aprés ce qui précéde. 

Nous cherchons, en fait, la variation de l’ampli- 
tude maximum (somme géométrique des compo- 
santes) en fonction de l’écart des fréquences w/w. 

En ordonnant suivant les puissances de wy/w 
léquation (4.6) on voit qu’en tenant compte des 
harmoniques 2 et 3, on ajoute a la composante 
en W,/w, de valeur absolue 2¢(a?+ b?)” une com- 
posante de valeur absolue (y/4)(a?+ 6?) et a la 
composante en (@,/w)? de valeur absolue 
(3x/4) (a? + b?)* une composante 


& £2 (a2 2h2\2 1 b2(202__ 2272 — % (q2 1 p2)*/2 
36 la" (a 3b?)? + b? (3a? — b?)?] 36 (4 4b)", 


2 2 2 2 
a) a) 3a w , a) a) 30a w : ; 
x= | 265 + 28 a—l) 38. 249 (at + b)| sin wt + E aot eet ria 5b (a2+ v)| cos wt + 
(4.6) 
a gah tel o x Wy 4 git: puitee 
peer tls — b?) sin 2wt + 2ab cos 2wt] 36 o2 a (a? — 3b) sin 3et + 6 (3a? — b?) cos 3t]. 
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En négligeant les phases relatives des diverses 
composantes, nous corrigerons donc l’équation 
(4.7) précédente, en remplacant les termes en ¢ 
et « par les nouvelles valeurs calculées, ce qui 
donnera une évaluation par excés de l’amplitude 
maximum (les composantes ajoutées étant sup- 
posées en phase avec la précédente). 

Remplagant comme précédemment (a?+ 62)” 
par x, il vient pour nouvelle équation de la figure 
de résonance: 


¥ 2 
(=+ Soe + day (+52) =1 (4.8) 


remplacant l’équation précédente (4.2). 

La différence ne porte que sur les valeurs 
numériques des coefficients « et y. 

Pour un calcul plus exact, il faudrait distinguer 
le voisinage de la résonance (1— w/w? négli- 
geable) et les zones éloignées de la résonance 


x 


(termes en « négligeable a cdté de 1—«4/w?). 


4.5, Autres propriétés de l’équation (4.1) 


On peut voir maintenant les effets physiques 
des coefficients en « et y de l’équation (4.1). 
L’introduction du terme d’élasticité 

linéaire «: 


non 


— diminue l’amplitude de résonance et aug- 
mente les efforts si «> 0. 


— augmente l’amplitude et diminue les efforts 
sia <0. 


L’introduction du terme y augmente l’effet de 
V’amortissement d’autant plus que l’amplitude 
est plus grande. Le maximum d’amplitude dimi- 
nue quand y croit. 

Enfin, calculons la courbe lieu des milieux des 
cordes verticales, qui remplace, dans le cas de 
cette figure de résonance, le diamétre conjugué 
de x= const. pour I’ellipse. 

L’équation de la courbe cherchée est: 


Tx ppl 
ya=x|l+ 9 x (e+ 5) | 


Au voisinage de l’origine, y équivaut a x(1—2é?), 
cette courbe est donc tangente au diamétre con- 
jugué de x=const. dans [ellipse représentant 
Voscillateur linéaire de méme amortissement, ou 
encore, la tangente a cette courbe a l’origine 
donne la valeur w, de la pulsation de résonance 
du systéme linéaire. 


5. Caleul d’un terme élastique non linéaire 
dans le cas de l’amortissement visqueux 
5.1. Expression de l’acuité de résonance 


Nous partirons de l’équation (4.1) sans terme 
en y pour laquelle la figure de résonance est 
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2 
(=+ ae + Axy 6? = (5.1) 
qui ordonnée en y devient: 
y+ Ixy (221-5 #) + 
(5.2) 


+ 42 (14 orale =—(). 


On calculera successivement le maximum 
d’amplitude (valeur double en y dans (4.9)) et la 
pulsation correspondante, puis les pulsations 
correspondant a une amplitude égale a la moitié 
du maximum, dont la différence est Aw. 

La discussion de |’équation (4.10), en fonction 
des paramétres, conduit au faisceau des courbes 
Fig. 7, limité a x> 0, présentant les caractéris- 
tiques suivantes: 


a>0 1 lobe fermé, 
a—=0 celobe est exactementl’el- 
lipse du systéme linéaire, 


1 lobe fermé plus une 


ele 2)2 
Saya ase Sectl branche infinie, 


9 
ae (lee 2 courbes se soudent 


ay en un point double, 
99 3 »\2 2 branches infinies sans 
Se Sage tae tangente paralléle a Oy. 


Les courbes correspondant 4 la méme valeur 
absolue de « se coupent sur y=x. II est facile de 
préciser le lieu des points 4 tangente paralléle a 
Oy, en annulant le discriminant de l’équation (5.2) 
puis en éliminant «: 


4e*x (y + e2x) —1= 0. (5.3) 


hyperbole centrée a l’origine, dont les asymptotes 
sont x=0 et y=— ex. 

En éliminant au contraire x*, on obtient la 
valeur de y/x correspondant au maximum: 


bh (2) =3 (l 36) 45 a—+5a] - (5.4) 


x Wo 


En reportant dans (4.10) la valeur (unique, réelle 
et positive) de l’amplitude x égale a la moitié du 
maximum, on obtient une équation du deuxiéme 
degré en y/x dont la résolution est longue mais 
sans difficulté. 

La différence relative des deux pulsations cor- 
respondant a l’amplitude moitié du maximum 
est alors: 

Aw 35 x 


reae /2 2 ae ae 5 
ae = 243 [b+ +5 mes | (5.5) 
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avec les approximations: « petit a cété de «?, et 
é=+ 1, permettant de remplacer les radicaux par 
des développements limités. 
Expérimentalement, on ne connaitra pas wp; 
pulsation du systéme linéaire non amorti, mais 
par exemple @,, pulsation du maximum d’ampli- 
tude ainsi qu’on a déja remarqué au paragraphe 
2.3. 
Avec les mémes approximations que plus haut, 
Péquation (4.12) conduit a: 
a5 a ep a 
oe hme it ea) e?(1—e?) 
et en faisant le rapport Aw/w, avec l’équation 
(5.5) on obtient: 
SeBlony/s [1+ 3et+3 3 ata | (5.6) 
on * e%(1—e?)|" * 
Cette équation dans laquelle 4w/w, est un résul- 
tat expérimental, permettra de calculer le terme 
délasticité non linéaire «, si l’amortissement ¢ 
est connu par ailleurs. Elle se raméne 4 Aw/w,= 
22/3 si x=0 et si « est petit. 
Expérimentalement, on connaitra, suivant 
lappareil de mesure utilisé, la pulsation du maxi- 
mum d’amplitude, ou de vitesse, ou d’accéléra- 
tion. 


5.2. Calcul numérique du terme élastique non 
linéaire 

A titre d’exemple, considérons le cas d’une 
barre d’acier de 400 x 18 x 6 mm? travaillant en 
flexion sur des appuis simples placés aux nceuds 
du premier mode de vibration, pour laquelle on a 
trouvé, autour d’une certaine amplitude (de 
Vordre de 300 um en son centre): 


— décrément logarithmique des _ oscillations 


libres: 6= 27ze— 0,24 - 10-3 
— acuité de résonance expérimentale: Aw/w,= 


0,38 - 10-2 


la premiére donnée conduit a: ¢=0,038- 10-3 
(le caleul de ¢ d’aprés la formule simplifiée con- 


duirait par la deuxiéme donnée a ¢ = ——-— 
Ot ireiO—*). 
On calcule ainsi « d’aprés l’équation (5.6) et on 
trouve «=12-10~°*. Ce calcul n’est donné qu’a 
titre indicatif. On ne peut évidemment mesurer 
que les coefficients du systéme oscillant formé par 
l'ensemble de la barre et des supports, avec tous 
les frottements et les pertes extérieures, et pour 
une amplitude donnée. Parmi ces pertes se trouve 
la dissipation magnétique dans |’électro-aimant 
d’excitation. Celui-ci était fermé sur la méme 
impédance au cours des deux expériences. Les 
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approximations faites devraient, d’autre part, 
étre revues de prés. I] serait enfin a discuter dans 
quelle mesure la relation = 2z¢ est valable pour 
un systéme non linéaire. 

Le calcul numérique fait ici a cependant le 
mérite de mettre en évidence influence sur la 
courbe de résonance, d’un terme non linéaire 
méme aussi faible que celui-ci. 


6. Représentation en coordonnées logarithmiques 
6.1. Principe de la transformation, cas de systéme 
linéaire 
Cette représentation, suggérée par M. J. Briz- 
LOUIN, consiste 4 tracer au lieu de l’ellipse de 
résonance précédente, la courbe transformée par 
é=logx 
n= logy. 
Dans cette représentation (Fig. 8) les droites 


x=(w/w,)2—=const. sont remplacées par des 
DA 0 Pp Pp 
droites a 45° 7=€-+ const. soit: 


(6.1) 


n= 2logw—2loga, + &. 


Fig. 8. Figure de résonance d’un oscillateur linéaire en 
coordonnées logarithmiques. 


Ces droites coupent l’axe 7 en des points que 
l’on peut graduer en @/ , lorigine correspondant 
a W=W, (pour w= nw, n=&+ 2logn). 

Les ellipses sont remplacées par une famille de 
courbes symétriques par rapport a 7=€, et ayant 
pour asymptotes 7—0 et €=0. 

Les maxima de x et y correspondent toujours 
(pour x et y positifs) aux maxima de € et 7,"et 
sont données par les tangentes paralléles aux axes. 

L’intérét de cette représentation est que la 
courbe expérimentale étant x«=f(w) on connait 


222 


y=(w/m>)?x au facteur wo” prés; donc, en co- 
ordonnées logarithmiques, la courbe (7, &) est 
définie 4 une translation prés. Connaissant sa 
symétrie, il est facile alors de la replacer dans ses 
axes. 


6.2. Systéme non linéaire 


En appliquant la méme transformation aux 
figures de résonance d’un systéme non linéaire, 
on obtiendra des courbes analogues aux précé- 
dentes, mais déformées et non plus symétriques. 

Soit, en effet, dans le plan x, y, l’équation (4.2): 


[ f(x) —y]?+ 4xy (e+ yx)?=1. 


Aprés transformation l’asymptote verticale de- 
vient: 7=— ©, &=log%», x, étant la valeur qui 
satisfait f(x))=1. L’asymptote verticale n’est 
donc plus = 0 sauf cas particulier, mais elle est 
connue expérimentalement, si la fonction f est 
connue. 

L’autre asymptote (x=0, [f(0)—y]?=1) sera 
comme précédemment =—o, 7=0 dans le 
cas le plus général ou f(0) est nul. 

La vitesse vibratoire maximum est toujours 
u=ox=yoo/o donc u? =xyw;. Les droites 
&+m=const. sont done des droites de vitesse 
constante, et la vitesse maximum est obtenue au 
point de tangente paralléle a 7=—€&. 

Pour un oscillateur non linéaire attaqué par 
une force donnée, on connait expérimentalement 
V’amplitude relative x (X, étant, en fait, trés peu 
différent de amplitude statique) et une varia- 
tion de wy) ne fait que translater la figure de 
résonance. 

Il est donc possible, a partir d’une courbe de 
résonance expérimentale de construire la figure 
de résonance correspondante en axes logarith- 
miques, puis, ayant suffisamment extrapolé vers 
les branches infinies, de placer approximative- 
ment les asymptotes, ce qui détermine, d’une 
part, la valeur de logx, tel que f(x))=1, d’autre 
part, l’axe O€ a partir duquel on peut graduer en 
@/@, les axes obliques a 45°. 

Pratiquement, il faut mesurer avec une préci- 
sion suffisante quelques amplitudes pour des fré- 
quences beaucoup plus élevées que celle de 
résonance. 


7. Autre représentation: le cercle de résonance 


7.1. Définition du cercle de résonance d’un oscilla- 
teur linéaire 


L’ellipse de résonance précédente ne permet 
guére une représentation graphique quand l’amor- 
tissement est trés faible, et surtout la zone entou- 
rant la résonance est trés mal représentée. 
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Il vient naturellement 4 l’esprit, d’abord de 
rapporter cette ellipse a ses axes principaux, puis 
d’effectuer une projection pour la ramener 4 un 
cercle. 

Soient OX, OY les nouveaux axes: (Fig. 9). 


L’équation (2.3) devient alors: 
2e7X 24+ 2Y2(1—e«?)—1=0 


que l’on raménera a un cercle, en multipliant les 
abscisses par (1—e?)”/e soit environ l/e si ¢ est 
petit: . 


2X242Y2=1/(1—e)~1. (7.1) 


La représentation des amplitudes est donnée par 
les droites paralléles a la tangente au cercle de 
pente (1—e?)”*/¢ ~1/e. 

Le rayon R du cercle mesure |’amplitude maxi- 
mum. I] peut étre gradué linéairement de 0 a 
100% en valeurs relatives, valeurs reportées sur 
des cercles concentriques (Fig. 9). La représenta- 
tion des pulsations est donnée par la pente 
d’un rayon vecteur, indiquée par la graduation 
de la tangente X= R au cercle. 


x 
fan] 
pa ~~ Resonance 
og damplitude 
este 
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Fig. 9. Cercle de résonance pour un systéme linéaire. e étant 
connu, recherche de la pulsation d’amplitude maxi- 
mum et construction d’un point quelconque. 


Le changement d’axes montre qu’on obtient, 
au moins pour les faibles amortissements, et au 
voisinage de la résonance, une graduation 
linéaire en pulsation de axe X=R, de la 
forme: 


¥=R(2—1) 0-8" ~ 2 (21) 


e € \w 
indiquant sur l’axe OX, m=, , et pour Y=R 
w/wmy=1+e(1—e?)* w1t+e. 


Ce diagramme sera done constant quel que 
soit e. 
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Les valeurs de w/q , dans le cas de résonances 
trés aigués, seront d’autant plus étalées que 
l’amortissement ¢ sera plus petit. Ayant défini e, 
on passe facilement 4 la courbe de résonance 
classique. 

Il est A remarquer que maintenant c’est |’am- 
plitude maximum qui sert d’unité, et non plus 
la déformation statique (amplitude pour w= 0). 

Si on voulait réaliser un abaque il semblerait 
d’ailleurs plus aisé de tracer d’avance des droites 
perpendiculaires 4 OX, et les rayons vecteurs cor- 
respondants, passant par les points d’intersection 
de ces droites et du cercle, et de faire tourner la 
droite portant les graduations de w/a , en fonc- 
tion de la valeur de «. 

La transformation inverse (courbe de réso- 
nance classique —cercle) est plus difficile. Elle ne 
présente d’ailleurs aucun intérét pratique quand 
le seul paramétre variable est e. 


7.2. Représentation dans les axes X Y d’un systéme 
non linéatre 

Au lieu de partir de l’équation (2.3) de lellipse, 
nous effectuerons les changements de variables 
précédents sur l’équation (4.8) ou pour simplifier 
sur (4.2) en sachant qu’on doit y remplacer en- 
suite « par 7x/9 et y par y/8. 

Afin de comparer a l’équation précédente du 
cercle, nous utiliserons les coordonnées polaires, 
avec X=ocos), Y=gsin6@ et en posant de plus 
€= cosq nous obtenons alors avec des transforma- 
tions trigonométriques simples: 


8yV20? . : 

2 0 2 — tag! 

207+ no sin? (y — 6) sin (p + 6) 
AxotsinOd . , sth 
ea OR ac 


On voit que la courbe correspondante sera trés 
voisine du cercle si y et ~ sont petits. 
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Il est en fait trés difficile de passer de la courbe 
de résonance expérimentale d’un systéme 4a cette 
représentation, car la graduation des pulsations 
n’y est connue qu’en fonction du paramétre ¢ 
supposé initialement inconnu. 


8. Conclusion 


Cette représentation de la résonance a |’avan- 
tage de fournir des courbes trés expressives du 
phénoméne. Elle permet de plus un calcul par 
approximation de certains oscillateurs dont les 
paramétres seraient connus. 

Enfin elle met en évidence le caractére non 
linéaire de certains oscillateurs expérimentaux. 

Une deuxiéme étape a franchir est d’évaluer 
les paramétres d’un oscillateur en partant de sa 
courbe expérimentale. L’étude de ce point est a 
peine ébauchée. I] nous a paru bon cependant 
d’exposer les principes de cette représentation, 
méme avant d’en avoir tiré les applications les 
plus souhaitables. 


Je remercie particuliérement M. H. Gonpvet, 
Directeur général du groupe de laboratoires de 
Bellevue du C.N.R.S., pour ses conseils précieux 
et ses encouragements, ainsi que M. J, BRILLOUIN, 
ingénieur-conseil, pour ses avis et ses suggestions, 
et spécialement pour la représentation logarith- 


mique du chapitre 6, qui lui est entiérement due. 
(Regu le 20 Aoat 1952.) 
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SCHALLDURCHGANG DURCH DEN MENSCHLICHEN KNOCHEN 


Von FRANZISKA SEIDL 


I. Physikalisches Institut der Universitat Wien 


Zusammenfassung 


Die beim Schalldurchgang durch Knochen auftretende Absorption und Reflexion wird mit 
einer Schlierenoptik sichtbar gemacht. Mit der Prismenmethode wird die Schallgeschwindigkeit 


im Knochen bestimmt. 


Summary 


The absorption and reflection of sound incident on human bones is made visible by means of 
Schlieren methods. Using a prism device, the velocity of sound in bone is measured. 


Sommaire 


On a rendu visible, au moyen d’un procédé optique de striation, l’absorption et la réflexion 
que subit le son quand il traverse les os. On a déterminé par la méthode du prisme la vitesse du 


son dans les os. 


Die Anwendung des Ultraschalls auf dem Ge- 
biete der Diagnostik pathologischer Veranderun- 
gen am Hirn setzt die Kenntnis der Vorgange 
beim Schalldurchgang durch Knochen voraus. 

Es wurden bereits von mehreren Forschern 
hinsichtlich der Reflexion, Absorption, Brechung 
und Fortpflanzungsgeschwindigkeit des Schalls 
im Knochen Angaben gemacht, die teils auf ex- 
perimentellem Wege, teils durch Abschatzungen 
auf Grund bekannter SchallgréBen gewonnen 
wurden. 

Ganz besonders méchte ich die Arbeit von 
W. Gtrrner, G. FrepLER und I. PArzoxp [1] 
hervorheben, in der alle Effekte aufgezeigt wer- 
den, mit denen bei der Durchschallung des 
menschlichen Schadels zu rechnen ist. Ferner 
wird in derselben, unter Beriicksichtigung der 
Gesetze der Strahlenoptik, der Strahlengang eines 
Ultraschallbiindels durch den Schadelknochen 
an Hand einer Skizze erldutert. 

Unabhangig von der Arbeit der obgenannten 
Verfasser wurden Untersuchungen an verschie- 
denen Knochen hinsichtlich des Schalldurch- 
ganges ausgefiihrt, wobei die Schlierenmethode 
angewandt wurde. Letztere findet man, infolge 
ihrer Unentbehrlichkeit fiir die Beobachtung der 
Ausbreitung der Schallwellen und der Struktur 
des Schallfeldes, so haufig in der Literatur ange- 
geben, daf} sich eine Beschreibung eriibrigt [2]. 
Doch sei bemerkt, daB auf die Einstellung des zu 
durchschallenden Praparates besondere Sorgfalt 
verwandt werden mu. Es kénnte z. B. der Fall 
eintreten, daB infolge einer Brechung die Schall- 
strahlen ungiinstig zur Richtung der Licht- 


strahlen geneigt sind und demzufolge entweder 
iiberhaupt nicht oder nur sehr schwach in Er- 
scheinung treten. Um dieser Unsicherheit auszu- 
weichen, wurden die Knochenpraparate um eine 
vertikale Achse drehbar montiert und auBerdem 
konnte die Stellung der Kiivette mittels Stell- 
schrauben verdndert werden. Beim Knochen- 
prisma war nur die Verdrehung um eine vertikale 
Achse von Bedeutung, da der schallabstrahlende 
Quarz auf der einen Prismenflache satt auflag. 

Abb. 1 zeigt den Durchtritt des Schallstrahlen- 
biindels durch die Scapula. Die Frequenz betrug 
800 kHz. Als Schallgenerator diente ein medizi- 
nisches Ultraschallgerat, auf dessen Schallkopf 
ein Tubus, der mit Wasser gefiillt und oben mit 
einer Gummikappe abgedeckt war, aufgesetzt 
wurde. Man erkennt, das das aus der Gummi- 
kappe austretende Schallstrahlenbiindel bereits 
auBerhalb des Interferenzfeldes vor dem Quarz 
liegt, und daB es infolge Reflexion der Schall- 
wellen am Knochen zur Ausbildung stehender 
Wellen kommt. Es werden somit im vorliegenden 
Fall die Schallwellen einerseits am Knochen 
reflektiert, andererseits erkennt man aber, dah 
ein betrachtlicher Teil der Schallenergie durch 
den Knochen hindurchgeht. 

Bei Schallabsorptionsmessungen, d.h. bei Mes- 
sung der Intensitat des Schallstrahlenbiindels vor 
und nach der Durchdringung des Materials, ver- 
meidet man das Zustandekommen stehender 
Wellen. Hier wurden sie mit Absicht mittels eines 
Reflektors erzeugt, um ein Kriterium dafiir zu 
haben, daB die Aufhellung im Dunkelfeld nicht 


von einer eventuell mangelhaften Ausbildung des- 
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selben herriihrt, sondern tatsachlich von dem 
austretenden Schallstrahlenbiindel. Man bemerkt, 
daB die Schallintensitat des Schallstrahlenbiindels 


dort abnimmt, wo die Knochendicke zunimmt. 


Es ist aber noch zu beriicksichtigen, da die 
Durchlassigkeit auch von dem Verhiltnis der 
Plattendicke zur halben Schallwellenlange ab- 
hangt, ferner von der Schwachung, gegeben durch 
die Streuung und wahre Absorption. Auch sei die 
Inhomogenitat des Knochenaufbaues hervor- 


gehoben. 
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nationalen UltraschallkongreB: in Rom (1950) 
aufmerksam gemacht. 

In vorliegender Arbeit wurde sie zur Bestim- 
mung der Fortpflanzungsgeschwindigkeit des 
Schalles im menschlichen Knochen angewendet. 
Aus einem mdglichst homogenen Knochenstiick 
wurde ein Prisma so herausgeschnitten, daB das 
Anschleifen der unter dem brechenden Prismen- 
winkel geneigten Flache an kompakter Stelle er- 
folgte. Wie man der Abb. 3 entnimmt, wurde ein 
kleiner Schwingquarz, dessen obere und untere 


Abb. 1. Abb. 2. Abb. 3. 
Durchschallung der Scapula. Durchschallung des Schallbrechung in einem 
Schadelknochens. Knochenprisma. 


Die Abb. 2 zeigt den Durchtritt des Schall- 
strahlenbiindels durch ein Knochenstiick der 
Calvaria, wobei die gleichen SchallgréBen hin- 
sichtlich Frequenz (800 kHz) und Intensitat 
(4 W/cm?) wie im vorherigen Versuch verwandt 
wurden. Bei diesem Praparat erfolgt die Schall- 
einstrahlung zufolge der gekriimmten Flache 
unter verschiedenen Winkeln. ErfahrungsgemaB 
ist die Durchlassigkeit dort am gréBten, wo die 
Einstrahlung in Richtung des Kriimmungsradius 
vorgenommen wird. Auch hier scheint sich dies 
auszuwirken. Allerdings ist auch eine geringe 
Dickenabnahme an der Durchtrittsstelle zu be- 
merken. Hinsichtlich der Inhomogenitat des 
Schadelknochens sei auf eine Abbildung in der 
oben zitierten Arbeit [1] aufmerksam gemacht. 

In Abb. 2 erkennt man wieder die Reflexion 
am Knochen. 

Die Strahlbrechung, die Abnahme der Schall- 
intensitat sowie die Fortpflanzungsgeschwin- 
digkeit des Schalls im Knochen — und damit bei 
bekannter Dichte den Schallwiderstand — kann 
man auf sehr einfache Weise mittels eines 
Prismas bestimmen, das aus dem hinsichtlich 
dieser GréBen zu untersuchenden Material gefer- 
tigt ist. 

Auf diese Methode wurde von J. SCHAFER [3] 
und spater von A. Barone [4] auf dem Inter- 


Flache vergoldet ist, mit einer Feder gegen die 
ebenfalls mit diinnster Goldschicht belegte Pris- 
menflache leicht angedriickt. Das Praparat 1aBt 
sich durch entsprechende, bereits erwahnte Vor- 
kehrungen sehr exakt einstellen, so daB gute 
Schlierenbilder erhalten werden, die eine Bestim- 
mung der Richtung der Strahlungscharakteristik 
zulassen. 

Die exakte Einstellung des Knochenprismas 
zum Lichtstrahl erkennt man an dem ausgebil- 
deten Beugungssaum der in der Abb. 3 als weibe 
Linie den Knochen, die Schrauben und die Feder 
umrandert. 

Das Schallstrahlenbiindel wird senkrecht zur 
Quarzoberflache abgestrahlt und trifft unter 
einem Einfallswinkel «x, der seiner GréBe nach 
dem brechenden Prismenwinkel gleich ist, die 
gegeniiberliegende geneigte Prismenflache. 

Mittels eines Zyklometers wurden der bre- 
chende Prismenwinkel als Einfallswinkel « und 
der Brechungswinkel /, der von der Richtung der 
aus dem Prisma austretenden Strahlungscharak- 
teristik und dem Einfallslot gebildet wird, ge- 
messen. Aus dem Snelliusschen Brechungsgesetz 
erhalt man fiir die Fortpflanzungsgeschwindigkeit 
der Longitudinalwelle im Prismenmaterial: 
sina 


ATS de 


226 


V bedeutet die Fortpflanzungsgeschwindigkeit 
des Schalls in der zur Beobachtung der Strah- 
lungscharakteristik gewahlten Flissigkeit, die im 
vorliegenden Falle CCl, war. Die Frequenz betrug 
3 MHz. 

Da infolge der Neigung der Prismenflachen 
die von den Schallstrahlen im Knochen zuriick- 
gelegte Weglange zunimmt, wird von einer be- 
stimmten Eintrittsstelle an, infolge Reflexion, 
Streuung und wahrer Absorption, die Schall- 
energie so geschwacht sein, da sie fiir die An- 
wendung der Schlierenmethode nicht mehr aus- 
reichend ist. 

Da aber hier nur auf eine Intensitatsabnahme 
hingewiesen werden soll und vorlaufig auf quanti- 
tative Angaben verzichtet wird, eriibrigt sich die 
Anwendung empfindlicherer MeBmethoden. AuBer 
den eben angefiihrten Griinden ist aber auch noch 
zu beriicksichtigen, dafs der Durchtritt der 
Schallwellen nur an denjenigen Stellen erfolgt, an 
denen die Dicke ein ganzzahliges Vielfaches der 
halben Schallwellenlange ist. 

Der relative Fehler 

l 


AV 
V; = cotg fh: AB 


betrug bei Ausmessung der Aufnahme Abb. 3 
mehr als 5%. Es wurde gefunden, dai} die Fort- 
pflanzungsgeschwindigkeit des Schalles an einer 
ziemlich kompakten, homogenen Stelle des 
menschlichen Knochens 3,8-10° cm/s betragt, 
wenn mit einer Schallfrequenz von 3 MHz ge- 
arbeitet wurde. Die Fortpflanzungsgeschwindig- 
keit des Schalles in Tetrachlorkohlenstoff wurde 
einer Tabelle entnommen. Eine exakte Messung 
derselben unter Beriicksichtigung der Tempera- 
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tur ware zwecks Erhéhung der MeBgenauigkeit 
hinsichtlich der Prismenmethode sehr erwiinscht 
gewesen, konnte aber nicht durchgefiihrt werden. 

Interessant ist, daf} das austretende, gebro- 
chene Strahlenbiindel von parallel verlaufenden 
Streifen durchzogen ist, was vermutlich auf die 
Ausbildung einer Kombinationswelle durch Re- 
flexion an den Kiivettenwanden zurickzufiihren 
ist. Solche Streifung zeigt auch die Strahlungs- 
charakteristik, die J. ScHArerR bei Anwendung 
der Keilmethode zur Bestimmung der Fortpflan- 
zungsgeschwindigkeit in Festkérpern erhalten 
hat. 

Es scheint empfehlenswerter zu sein, nicht das 
zentrale Bild, sondern samtliche Beugungsbilder 
bei der Schlierenmethode abzublenden, wenn 
man die Aufnahme ausmessen will [3]. Aus den 
Abbildungen in der Arbeit von A. BARONE [4] 
geht hervor, da es vorteilhaft ist mit hohen 
Frequenzen zu arbeiten, da sie eine gute Biinde- 
lung der Schallstrahlen bewirken. 

Die hier mitgeteilten Versuchsergebnisse zei- 
gen, das Schallstrahlen an Knochen reflektiert 
und beim Durchgang gebrochen werden. Dieses 
Verhalten mu bei Anwendung des Ultraschalls 
zur Abbildung am menschlichen Schadel zu Bild- 
verzerrungen fiihren, was bei der Dussixschen 
Methode zu beriicksichtigen ist. 

(Eingegangen am 17, Januar 1953.) 
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THE ATTENUATION OF SOUND PROPAGATED 
OVER THE GROUND 


By J. D. HAYHURST 


Ministry of Supply, London 


Summary 


Theoretical calculations by KNUDSEN, supported by the results of laboratory experiments, 
have established the attenuation of sound in still air. Little work has hitherto been done out-of- 
doors because of the difficulty of making allowances for the effect of wind on the attenuation. In 
the course of an aircraft noise abatement investigation made recently by the Ministry of Civil 
Aviation, the effect of wind on the propagation of sound over a dry concrete runway has been 
explored up to a distance of about half a mile from a source of sound. It was found that the only 
significant parameter was the component of wind in the direction of propagation and the values 
that emerged showed that wind effects cannot be neglected in any acoustic work made out-of- 
doors. By interpolation the attenuations corresponding to a zero wind component in the direction 
of propagation were derived and were found to be substantially greater than those previously 
determined for still air. The attenuations were however statistically independent of the absolute 
wind speed. Repetition of the investigation but over grassed areas produced no reliable results. 


Sommaire 


Au cours d’une étude sur la réduction du bruit d’avion, on a déterminé, jusqu’a une distance 
d’environ 800 m d’une source sonore, l’influence du vent sur la propagation du son au-dessus 
d’une piste séche en béton: seule la composante du vent dans la direction de la propagation a une 
influence. On a déduit par extrapolation de ces résultats les affaiblissements relatifs 4 une com- 
posante nulle du vent dans la direction de propagation: ils sont nettement plus grands que ceux 
trouvés précédemment en air calme. Ils sont cependant statistiquement indépendants de la vi- 
tesse absolue du vent. En répétant les essais au-dessus de zones herbeuses, on n’a pas obtenu de 
résultats certains. 


Zusammenfassung 


Die Berechnungen von KNUDSEN iber die Schallabsorption in ruhender Luft sind bisher fast nur 
in Laboratoriumsversuchen bestatigt worden, weil der Einflu8 des Windes bei Messungen im 
Freien nur schwer zu eliminieren ist. Bei emer Messung der Schallausbreitung iiber einer trockenen 
Betonrollbahn wurde festgestellt, daB der einzige wesentliche Parameter die Komponente des 
Windes in der Ausbreitungsrichtung des Schalles war. Bei der Extrapolation der Schallabsorption 
auf verschwindende Windkomponente in Ausbreitungsrichtung wurde ein erheblich héherer Wert 
gefunden, als friiher fiir ruhende Luft bestimmt worden war. Die Absorption war jedoch stati- 
stisch unabhangig von der absoluten Windgeschwindigkeit. Messungen iiber grasbewachsenen 
Flachen ergaben keine zuverlassigen Resultate. 


1. Introduction 


During 1951 the Ministry of Civil Aviation 
made an investigation to explore the noise abate- 
ment pruperties of bays to contain aircraft, the 
engines of which were being tested. The investiga- 
tion consisted, in essence, of running an aero- 
engine behind an experimental wall of asbestos 
cement sheeting progressively extended in area 
and measuring the sound pressure levels in front 
of the wall. The reduction of sound pressure levels 
below those when there was no wall interposed 
gave, at various distances from the wall, the 
screening power associated with a particular size 
of wall. The results of the investigation, which 
has been reported elsewhere [1], led the Ministry 
of Civil Aviation to construct a prototype run- 
ning-up pen at London Airport. 


Before the investigation proper commenced, 
attention was given to the variables which would 
affect a comparison of the screened and unscreened 
sound pressure levels and a series of measurements 
was made to find the way in which the sound 
pressure levels varied with the distance away 
from the source of sound and how this variation 
depended on atmospheric conditions, in particular 
the wind. 

This paper describes the measurements that 
were made and the conclusions that were drawn 
from them. 


2. The test programme 


The tests were planned against the background 
of the principal investigation, in which sound 
pressure levels were to be measured over a range 
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of distances up to about 3000 feet from the ex- 
perimental wall being built on a disused runway 
at London Airport.* In order to eliminate the 
effects of varying terrain the measurements were 
to be confined, in the main, to the line of the run- 
way which was of level concrete in open country 
without any buildings or trees. The sound output 
of the aero-engine was found to be consistent 
within the errors of measurement, so that at a 
distance from it on a given day (and the principal 
investigation was expected to take three months) 
the sound pressure levels would be dependent 
only on the atmospheric conditions. 

The variables in the atmospheric conditions 
which would affect the propagation were taken 
as humidity, temperature, wind velocity and 
direction, together with temperature gradient 
and wind gradient. For these last two, no means 
were readily available for their measurement and 
their effects were therefore allowed to fall within 
the errors of the experiment. Corrections for 
temperature and humidity were obtainable from 
published work of KNuDSEN [2]. There was some 
doubt, however, whether corrections for wind 
would be necessary within the limited range of 
measurements of sound pressure levels, since the 
levels would not be expected to be measurably 
altered either by an apparent linear translation 
of the source by the wind or by a change of path 
distance by refraction by wind gradients. It was 
also accepted that the propagation would be af- 
fected by turbulence in the atmosphere and by 
scattering (cf. [3]). In view of these doubts, it was 
decided to make a series of measurements in dif- 
ferent atmospheric conditions and see whether 
results could be obtained which would permit 
empirical corrections to the applied in the prin- 
cipal investigation. 

The source of sound was provided by a Bristol 
Hercules 630 14-cylinder aero-engine driving 
through a reduction gear of 0.444:1 a four- 
bladed wooden propeller 13.25 feet in diameter. 
The engine, mounted in the airframe of a Vickers 
Viking aeroplane, was operated at the constant 
setting of 30 inches of mercury manifold pressure 
and 1900 r.p.m., and at this power setting de- 
veloped about 750 b.h.p. 

The sound pressure levels were measured on a 
Standard Telephones and Cables noise meter with 
electrical filters inserted to analyse the levels in 


* Conversion table of units used to metric units: 


1 foot = 0.305 m, 
1 b.h.p. = 745.7 W, 
1 knot = 0.515 m/s, 
1 dB/1000 ft = 3.28 dB/km, 
1 dB/1000 ft/knot = 6.37 dB/km per m/s. 
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the eight octaves between 37.5 and 9600 c/s. The 
microphone stood, except for one test, on a tripod 
about 4 feet above the ground and was surrounded 
by a wire cage covered with fine muslin to prevent 
the wind affecting its response. 

Meteorological data for the investigation were 
taken from the London Airport Meteorological 
Office which was situated about 1 mile from the 
test site. 


3. The measurements made 


Nineteen series of measurements were made of 
the sound pressure level per octave in each of the 
eight octaves, the measurements being continued 
up to the distance at which the level fell below the 
range of the equipment or was indistinguishable 
from the local noise background. As would be 
expected, the higher frequencies were attenuated 
the most and the measurements in these frequen- 
cies were necessarily limited. Observations were 
made at intervals of 100 feet up to a distance of 
1000 feet from the aeroplane, thereafter at 250 feet 
intervals up to 2000 feet and finally at 2500 feet. 
The effect of the greater attenuation at the higher 
frequencies was such that, on a calm day, measure- 
ments in the eighth octave could be made only at 
the first five observation points, whereas readings 
over the full distance were possible in the first four 
octaves. This limitation in the higher octaves 
was exaggerated by the lower content of the 
higher frequencies in the sound emission of the 
source; the sound spectrum 100 feet in front of 
the aeroplane is shown in Fig. 1. 


110dB ai atl aie 7 


S 100 


8075 ve 150 300 600 1200 2400 4800 
-5 -150 -300 -600 ~-1200 ~-2400 -4800 -9%00ck 
—— Frequency octaves 


Fig. 1. Sound spectrum of Vickers Viking 100 ft. ahead. 


4. Analysis of the results 


In selecting a method of analysing the results, 
account had to be taken of the order of the experi- 
mental errors inherent in the tests, and of the 
immediate purpose of the tests. Measurement of 
sound pressure levels in the open air is far from 
being a precise art, especially as the levels tend to 
fluctuate because of the wind. Furthermore, only 
some of the atmospheric conditions were known 
and so the results contained unknown factors. 
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However, the immediate application of the results 
was to the principal investigation and the aim 
therefore was to deduce corrections which could 
be used with some degree of confidence in that 
investigation. 

It was therefore assumed that the sound pres- 
sure levels decreased with distance, in part by the 
inverse square law and in part linearly with dis- 
tance. In the second term, the linear variation 
would depend on an attenuation coefficient de- 
termined by the terrain (in this case level con- 
crete) and by the atmospheric conditions. 

Thus if L, and L, are the sound pressure levels 
at distances S, and S, from a source 


L, = L,—20log = 


2 


M(S;—Sz) 
where M is the attenuation coefficient. 

From a set of measurements of I and S, it is 
possible to calculate a curve of the form repre- 
sented by this equation which gives the best 
statistical fit to the observations. 

This was done for each of the eight octaves for 
each series of measurements, and the correspond- 
ing values of the attenuation coefficients cal- 
culated. These values are given in Table I. 

The attenuation coefficients in each octave 
were then corrected to a relative humidity of 
50%, a value which was arbitrarily chosen, but 
the variation of temperature over the test period 
suggested that no correction for temperature was 
required. In fact, corrections for humidity might 
also have been neglected. 

The corrected attenuation coefficients were 
next examined in relation to the wind conditions 
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obtaining during each series of measurements. 
Using the fourth octave as an example, it was 
found that there was no significant connection 
between the attenuation coefficients and the wind 
velocity ignoring its direction, i.e. the scalar 
wind. There was, however, a significant connec- 
tion between the coefficients and the component 
of wind velocity in the line of observations, i. e. 
the vector wind. Fig. 2 is the scatter diagram for 
this fourth octave with the sign convention that 
a positive vector wind implies a wind component 
towards the source of sound. The regression line 
of attenuation coefficient upon vector wind has 
been included and the corresponding correlation 
coefficient is significant. 

This procedure was applied to the eight octaves 
and it was found that in five octaves the corre- 
lation coefficient was significant, but in the first, 


30 
dB/1000 ft 
20 


10 


—— Attenuation 


=f 
270 -10 0 
— Vector wind 


with 


10 knots 20 


Fig. 2. Variation of attenuation vector wind 


300—600 c/s octave. 


Table I. Calculated attenuations over concrete 


Attenuation in dB/1000 ft 


Relative | Scalar 
Set mgd Frequency octaves in c/s 
37.5 —75| 75 —150 | 150 —300] 300 —600 | 600 —1200 | 1200 —2400 | 2400 —4800 | 4800 —9600 
=15 == 1 — 1 0.5, + 1.5 + 7 

1 2.4 el 4.3 1.8 — 2.2 4.5 + 3.1 +32.9 
2 TON 2:5, 4-57 +. 9.9 +20.7 9.3 20.6 +52.3 
3 ae ais 3.2 +25.6 + 23.6 9.4 42.8 = 

4 8.5] + 2.0 +0.9 + 0.8 + 0.4 0.1 Dal +21.2 
5 0 — 0.6 | +4.1 + 7.3 9.6 9.0 8.5 - 33.4 
6 8.9 ahs) +1.1 ley 1-9 3.8 3.5 +38.9 
7 4.6] — 0.8) +0.8 + 1.2 4.2 3.8 6.1 +37.0 
8 4.4 7.4 3.6 4.6 5.4 4.4 till +19.4 
9 4.4. 13.4 + 1.2 +10.1 5.1 5.2 8.2 +33.0 
10 0 0 +4.1 + 1.6 2.4 6.3 - 3.8 + 26.0 
11 1.7] — 2.0 | —9.0 — 7.0 4.0 12.0 +23.0 17.0 
12 3.2] + 7.8 | +6.5 +10.5 3.0 15.6 + 7.9 +23.4 
13 3.2] +15.6 + 6.2 + 5.9 Le 17.1 +20.7 +40.4 
14 1.4 0 —4.0 + 7.0 3.0 4.0 4.0 +22.0 
15 7.0 1.9 | +7.9 +10.7 3.0 15.2 19.1 +38.0 
16 Bal 0.2 | +1.4 + 5.7 5.8 7.0 + 18.4 -48.4 
17 7.0 2.5 5.5 +12.6 3.4 + 24.9 +19.9 -41.0 
18 14.1 a5 2.6 2.1 1.8 — 1.2 + 0.1 +22.4 
19 5.0 1.9 | +5.0 + 4.2 0.5 — 0.1 — 1.3 —38.2 
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Table II. Effect of wind on attenuation over concrete 


Frequency octaves 


Attenuation in \ 
zero vector wind 


75 —150 | 150 —300 
+1.3 +5.0 


=N5 
+1.3 


Change of attenua- 
tion with wind 


+0.1 


Correlation coeffi- 
cient 


second, and eigth octaves it fell below any accept- 
able significance level. In addition, the values of 
the attenuation coefficients corresponding to zero 
vector wind were extracted and these are given in 
Table II which also includes the slopes of the 
regression lines, which give the rate of change of 
attenuation coefficient with vector wind. Having 
regard to the difficulties that have already been 
mentioned of measuring and analysing sound 
pressure levels out-of-doors, it was pleasing to 
find a significant correlation in five out of eight 
octaves. Of the remaining three octaves, the wind 
effects are apparently small compared with the 
experimental errors in the two lowest octaves, 
and the few sound pressure level measurements 
that could be made in the top octave did not lend 
themselves to an accurate determination of the 
corresponding attenuation coefficients. 

The attenuation coefficients corresponding to 
zero vector wind are given in Table II and are 
plotted against frequency in Fig. 3, from which 
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Fig. 3. Attenuation over concrete in zero vector wind. 


Ss -e 2.dB/1000 ft-knot 
Ss. 


2: 

2s 

= 8 

@< 

250 

cp foe Mee: 150 300 600 1200 2400 4800 
-75 ~-150 -300 -600 -1200 -2400 -4800 -9600cIs 


Frequency octaves 


Fig. 4. Change with vector wind of attenuation over con- 
crete. 


+9.2 


+13.5 +21.0 +28.6 | dB/1000 ft 
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it will be seen that the non-significant coefficients 
lie substantially on a curve passing through the 
significant coefficients. 

The rates of change of attenuation with vector 
wind are also given in Table II and are plotted 
against frequency in Fig. 4. In this case the non- 
significant values do not lie on the straight line 
which passes through the significant values. It 
should be noted that there is no a priori reason 
why the attenuation coefficients should depend 
linearly on vector wind, nor why this variation 
itself should depend logarithmically on frequency. 
In the context of this present investigation, how- 
ever, assumptions of linearity may be justified. 


5. Consideration of results 


These results show that over the range of dis- 
tances explored a significant association can be 
established between the attenuation coefficient 
and the vector wind. It does not, of course, follow 
that it is the vector wind that itself is causing the 
changes of attenuation since if something else is, 
then the vector wind may be correlated with that. 

The attenuation coefficients at zero vector wind 
differ from those established by previous workers. 
In Fig. 5 have been plotted the results of the 
present investigation (from Fig. 3), the theoretical 
results of KNUDSEN [4] (which are confirmed by 
laboratory measurements), and the experimental 
results of Eyrine [5]. It will be seen that the 
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present investigation gives higher attenuations 
in the lower frequencies although the results are 
more comparable in the higher frequencies. It is 


true that the present results were obtained over | 


the ground when some terrain absorption might 
be expected, although the absorption charac- 
teristics of dry concrete are unlikely to account 
for the differences. Another possible explanation 
is that the experimental results out-of-doors may 
include effects of atmospheric turbulence and it 
is unfortunate that no assessment could be made 
of the gradients of temperature and wind that 
occurred while measurements were being made. 

The variation of attenuation coefficient with 
vector wind is not inconsistent with everyday 
experience but is much larger than would be 
expected from a consideration of the theory of 
sound propagation, since the changes of attenua- 
tion that have been obtained correspond to an 
extremely large change of track distance. Conse- 
quently, it is not easy to account for the wind 
effect either in terms of a virtual transposition of 
the source or in terms of an extended path caused 
by refraction. Again, turbulence offers a possible 
solution, but the effect would have to be a func- 
tion of the turbulence along the line of propaga- 
tion only, since there was no correlation between 
attenuation and the absolute wind. 

The results, then, are intriguing. It has not yet 
been possible to make an exhaustive enquiry as to 
their physical explanations but such an enquiry 
would appear to be worth while. 

One feature, however, emerges very clearly. 
In making acoustic measurements out-of-doors, 
it is not sufficient to describe the wind conditions 
in vague terms such as “‘slight breeze”’. The effects 
of wind are such that experimental results may 
be distorted by the presence of any wind at all 
and grossly misleading conclusions may be reached 
if adequate corrections are not made for wind. 


6. The effect of height on attenuation 


As has previously been mentioned, the normal 
procedure throughout the investigation was to 
have the microphone set up on a tripod about 
4 feet above the ground. When the opportunity 
occurred two series of measurements were made 
with the microphone placed at different heights. 
In one series, when there was.a wind component 
following the sound, the microphone was sus- 
pended from the platform of a travelling servic- 
ing trolley. In the other, when there was a wind 
component opposing the sound, the microphone 
was suspended from the hook of a crane. In each 
case, the sound pressures were measured with the 


microphone at 0, 10, 20 and 30 feet above the 
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ground and over the usual range of distances from 
the aeroplane. 

Measuring downwind with a vector wind of 
5 knots, it was not possible to pick out any varia- 
tion of attenuation with height. Upwind, with a 
vector wind of 7 knots, however, there was a pro- 
gressive decrease of attenuation with height. The 
results obtained are shown in Fig. 6. It is, of 
course, not possible from this single test to make 
any generalisations but the result is of sufficient 
interest to warrant consideration of further 
experimental work. 
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Fig. 6. Change of attenuation with height. 


e—-—e 0 feet above the ground, 
x —+—-:—x 10 feet above the ground, 
gw ———— @ 20 feet above the ground, 
+r++++++ 30 feet above the ground. 


7. The effect of terrain on the propagation 
of sound 


In exploring the effects of wind on propagation, 
all the measurements were made with a con- 
tinuously paved surface between the test aero- 
plane and the points of observation. The oppor- 
tunity was taken of making a series of measure- 
ments with grassed areas interposed. The same 
technique as before was adopted, but the attenua- 
tions obtained (given in Table III) were now cor- 
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Table III. Calculated attenuations over grass 
Tempe- | Relative | Scalar | Vector Attenuation in dB/1000 ft 
Set| rature | humidity] wind wind Frequency octaves in c/s 
37.5 —75 | 75 —150 | 150 —300]| 300 —600 | 600 —1200 |} 1200 —2400 | 2400 —4800 | 4800 —9600 
1 10} 3.0 + 3.0 | + 6.2 +11.1 +23.4 +26.6 + 63.0 
2 2.04) +125 Ubi |) Se A0st + 5.1 Te + 2.9 +16.0 
3 0.7 1.4 0.8 4.9 ar +20.7 +18.9 +62.4 
4 20) dal 5.4] + 3.6 + 8.5 +17.9 +23.4 +40.0 
5 3.6] +2.0 ee ee eal + 7.4 + 9.6 + 28.0 +38.0 
6 1.1} +5.5 6.3 8.3 + 4.5 +15.3 +20.2 +11.6 
7 et | me hale Ce ap te? + 4.0 +—6.7 +36.4 
8 0.4) +0.5 1.8} + 5.2 + 7.2 +15.2 +26.0 +26.0 
9 9.4) +2.5 —- 1.6] + 3:2 +10.0 te Coli +10.8 +11.0 
10 0.5] +4.6 3.6 8.3 +13.3 + 14.2 +12.0 +12.0 
ll 10.2] +4.5 + 2.7 - 6.6 +18.1 + 10.6 + 4.1 +20.2 
12 14.0} +4.9 +12.7 LOSS + 24.6 + 26.6 +57.0 + 88.0 


Mean difference of attenuation under 


same atmospheric conditions 


To obtain 95% confidence limits add 


rected for the vector wind appropriate to each set 
of measurements and the differences of attenua- 
tion over concrete and over grass computed for 
the same standard atmospheric conditions. The 
results are given in Table IV, the mean difference 
being plotted in Fig. 7. 

These results are not convincing. It was not 
possible to find continuously grassed surfaces and 
the lines of observation were crossed by concrete 
runways and taxi-ways, the intersections being 
not really negligible compared with the ranges 
of distances covered. Nor were the grassed areas 
themselves uniform, varying from newly sown 
areas, established areas recently mown, and areas 
where the grass was about one foot high. Further- 
more, it is possible that the wind effect on pro- 
pagation over grass will differ from that on pro- 
pagation over concrete since the wind will change 
the profile of the grass and may affect the absor- 
bent properties of the terrain, thus invalidating 
the wind corrections. These shortcomings are 
reflected in the ranges between the confidence 
levels shown in Table IV so that the mean differ- 
ences have little significance. 
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Fig. 7. Mean differences between attenuations measured 
over grass and over concrete. 
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In these circumstances, the only conclusion 
that can be drawn is that the attenuation is 
greater over grass than over concrete, and the 
difference increases as the frequency of the sound 
increases. These general conclusions were al- 
ready known [6]. 


8. Conclusions 


Sound pressure levels have been found to fall 
off in the open air at a greater rate with distance 
than would be derived from a combination of the 
inverse square law and of atmospheric absorption 
determined under laboratory conditions. 

Corrections for the effect of wind on sound 
propagation have been established which can be 
applied with some confidence up to distances of 
about one mile from the source. Effects of wind 
cannot be neglected in acoustic measurements 
made out-of-doors. 
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ON EDGE TONES AND ASSOCIATED PHENOMENA 


by ALAN POWELL 


Department of Aeronautical Engineering, University of Southampton 


Summary 


The edge tone phenomena occurring at high jet speeds have been photographed by the Schlieren 
method and the sound waves made clearly visible. They are of opposite phase on each side of the 
stream, and microphone traverses have shown this to be the case at low speeds also. This observa- 
tion makes possible an explanation of the mechanism whereby the sound waves give rise to 
stream disturbances at the jet orifice. Two basic relationships are given, one connecting the relative 
phases of the motion and the other dealing with the effective efficiencies of each of the component 
processes. Although an exact theoretical solution is not yet possible, certain assumptions with 
regard to disturbance velocity result in good agreement with experiment, after invoking the 
condition that the jet shall operate near the region of maximum sensitivity. It is pointed out that 
the lower frequency limit of the sensitivity arises from an insufficiently large rate of amplification 
of stream disturbance, whereas the upper limit is the result of a trend towards stability. The 
general theory is capable of explaining the jumps of frequency, associated hysteresis effects and 
other characteristics. The observed susceptibility to local reflectors or resonators is discussed, and 
it is suggested that an extension to other phenomena, for example the organ pipe, could be made. 


Sommaire 


On a photographié par la méthode des stries les phénoménes acoustiques produits par une aréte 
placée dans un jet de grande vitesse et on a rendu nettement visibles les ondes sonores produites. 
Les phases de ces ondes de chaque cété du courant sont en opposition et une exploration faite en 
déplagant un microphone a montré qu’il en est de méme aux faibles vitesses du jet. Ce résultat 
permet d’expliquer le mécanisme suivant lequel les ondes sonores produisent des perturbations 
du courant a l’orifice. On donne deux relations fonda mentales, l’une entre les phases relatives du 
mouvement et l’autre concernant les rendements effectifs de chaque composante du phénoméne. 

La théorie générale permet d’expliquer les discontinuités de fréquence, les effets hystérétiques 
associés, etc... On étudie l’influence observée de réflecteurs ou de résonateurs locaux, et on suggére 
d’étendre cette étude a d’autres cas (tuyaux d’orgue). 


Zusammenfassung 


Die bei hohen Strahlgeschwindigkeiten auftretenden ,,Schneidenténe“ wurden mit einer 
Schlierenmethode untersucht. Die Schallwellen haben auf den beiden Seiten des Strahles verschie- 
dene Phasen, auch bei kléiner Strahlgeschwindigkeit. Diese Beobachtung ermoglicht eine Erkla- 
rung der Entstehung von Strahlstérungen an der Diisendffnung durch die Schallwellen. Obwohl 
eine exakte theoretische Lésung noch nicht méglich ist, ergab sich mit der Annahme einer ,,St6- 
rungsgeschwindigkeit gute Ubereinstimmung mit dem Experiment unter der Voraussetzung, daB 
die Strémung sich nahe dem Bereich héchster Empfindlichkeit befindet. Es wird gezeigt, da} die 
untere Frequenzgrenze dieser Stérungsempfindlichkeit durch ungeniigende Verstirkung der 
Strahlst6rung und die obere durch eine Neigung zur Stabilitét bedingt ist. Mit der allgemeinen 
Theorie lassen sich die Frequenzspriinge, damit verbundene Hystereseerscheinungen und andere 
charakteristische Vorginge erklaren. Weiterhin wird die beobachtete Empfindlichkeit gegen 
Reflektoren und Resonatoren diskutiert und eine Anwendung auf andere ahnliche Erscheinungen, 
wie z. B. bei der Orgelpfeife, vorgeschlagen. 


1. Introduction 


“The problem of edge tones is one which con- 
tinues to form a battleground for rival theories, 
though a complete solution seems as far off as 
ever.” 

Such were the opening words of a paper by 
LentHan and Ricwarpson in 1940 [1] and al- 
though some further work has since been done, 
the situation has not developed to any great 
extent. The problem presented by edge tones lies 
in explaining the mechanism whereby the placing 


of a wedge symmetrically in the stream issuing 
from a narrow slit gives rise to a characteristic 
note and, it has been observed, an associated 
vortex pattern in the flow. The subject is in fact 
remarkable for the diversity of the nature of the 
theories suggested in the past, which Brown [2| 
discusses in some detail. At present they can be 
conveniently grouped into three main schemes. 


These are, briefly: 


(i) The basic mechanism is essentially hydro- 
dynamic and fairly localised in the region of the 
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vortex formations around the wedge. A deflection 
of the stream to one side gives rise to a pressure 
which is supposed to disturb the approaching 
stream so as to cause it to move to the opposite 
side. Brown [2] is largely responsible for this 
theory, which has been amplified somewhat by 
Jones [3]. However, the scheme is somewhat 
lacking in concrete evidence, and an attempt has 
yet to be made at the admittedly difficult task of 
placing the theory on anything like a quantitative 
basis. 

(ii) The vortices forming in the jet stream 
(which are supposed to form a KARMAN street) 
are supplemented by others springing from the 
edge of the wedge, so giving rise to a pair of 
KARMAN streets, one on each side of the wedge. 
This gives rise to a pendulation of the jet, so 
causing the edge tones. The existence of the 
secondary vortices springing from the edge have 
also been observed by Brown [4], but it seems 
that they appear only under certain conditions. 
Ricwarpson [5], who suggested this “double 
street” theory, found the vortex velocities to be 
near those estimated from RAYLEIGH’s result [6], 
rather than the approximately constant value 
one would expect for a KARMAN street, and the 
initiation of these vortices has still to be ex- 
plained. Thus this theory too, has its difficulties, 
and more work is necessary to substantiate or 
refute it. 

(iii) The sound wave produced by the vortices 
striking the edge creates embryo vortices at the 
orifice. Several attempts have been made to de- 
velop this theory, but without great success (see 
Brown [2] for a note on these). While noting that 
the whole of the edge tone phenomenon takes 
place within the range of sensitivity (to sound) of 
jets, BRown only makes use of this fact during 
the starting process of the theory he suggests. 
RicwHarpson [5] has reported an experiment 
which bears heavily against this class of mecha- 
nism. In this, two jets were set fairly close to- 
gether, only one of them being provided with an 
edge. While the latter produced an edge tone in 
the usual way, it executed no periodic influence 
on the one having no edge. 

The fact that the edge tone phenomenon takes 
place within the range over which jets are found 
to be sensitive to sound, despite several “‘jumps” 
in frequency of the edge tone, cannot be dis- 
missed lightly, since the edge tone is itself a 
source of sound located close to the orifice. It is 
the writer’s aim in this present paper to show 
that there is considerable evidence favouring a 
mechanism based upon the last of the groups 
above, that is depending upon and closely as- 
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sociated with the so-called sensitivity of jets. 
Brown [4] incidentally, considers the vortex 
motions of sensitive jets and edge tones to be 
similar. Perhaps it is appropriate to mention here 
that it will become evident that RICHARDSON’s 
experiment, mentioned above, is not as decisive 
as appears at first sight. 


2. Experimental findings 


During the execution of a research involving 
the optical study of highspeed jet flows (POWELL 
[7]), a phenomenon displaying certain character- 
istics found in sensitive jets and edge tones was 
revealed, although the Reynolds number was far 
higher than customary in these studies. 

The method used was a two-mirror Tépler 
Schlieren system, some slight refinements of ar- 
rangement and technique enabling sufficiently 
intense sound waves to be photographed. The 
basic scheme has been described in many publi- 
cations (see for example BARNES and BELLINGER 
[8] where a bibliography is given), a spark light 
source having an effective duration of about a 
microsecond being used. It is only necessary to 
comment here upon the salient characteristic of 
the system: it responds to the total density gra- 
dient in a direction perpendicular to a knife edge 
(the attitude of which is adjustable) taken along 
traces normal to the field shown on the photo- 
graphic plate. Density gradients in one direction 
are depicted by an increase in illumination, and 
in the opposite sense by a corresponding decrease, 
the system being quite insensitive to gradients in 
the direction of the knife edge. 

The streams with which we were concerned 
issued from a small rectangular orifice, at a speed 
equalling the local velocity of sound and with 
a static pressure at the exit exceeding that of the 
local ambient atmosphere, so giving rise to the 
characteristic regular or cellular pattern associat- 
ed with such jets (Pack [9], Harrmann and 
Lazarus [10]). It was found that disturbances 
in the jet, having an asymmetric (or alternate) 
pattern, gave rise to sound waves of definite fre- 
quency and of opposite phase on each side of the 
jet. A theory of the mechanism responsible for 
stabilizing the motion at a frequency uniquely 
associated with a flow parameter has been sug- 
gested (Powe. [11]), in satisfying agreement 
with experiment. The sound produced by inter- 
action of the stream disturbances with the regu- 
larly spaced shock waves of the jet is supposed 
to result in the creation of boundary layer dis- 
turbances at the jet exit. These are assumed to 
become amplified as they pass along the stream 
and so in turn are responsible for a sound output. 
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With regard to “sensitive” jets and the tones 
emitted by choked jets, an important point arises. 
The fact that the sound waves are in antiphase 


across the jet makes “‘sensitivity”’ possible, since: 


pressure gradients across the stream are thereby 
generated. This contrasts with the usual pheno- 
menon when an external source of sound is em- 
ployed, there then being no response when the 
sound is directed along the stream, as the forces 
tending to produce transverse deformations of 
the stream are then a minimum. 

Also, the sound effectively originates from a 
particular region of the jet, as in edge tones where 
it is generally understood that the sound ema- 
nates from the region of the edge. Thus in the one 
case the frequency is governed by the shock po- 
sitions and the state of the boundary layer, 
whereas in the other the edge does this duty. 

It thus becomes very interesting to determine 
whether the edge tone phenomenon can be set up 
in such a jet. To investigate this a bi-convex aero- 
foil of 2.5 em chord was introduced into the 
stream (emerging from an exit 3 mm wide) so 
as to form the edge. 


Fig. 1. Edge tone at high speed with Schlieren knife-edge 
parallel to jet. 


Two typical photographs are shown in Fig. 1 
and 2, the usual precautions against vibration of 
the knife-edge and the like being taken. It is at 
once clear that motion characteristic of edge tones 
is present, although the disturbances do not de- 
velop into the discrete vortices as may be some- 
times observed in the low speed case. In the first 
photograph, a knife edge parallel to the stream 
was used, and bearing in mind the fact that as the 
sound waves are proceeding in opposite directions 
on each side of the stream, it is clear that they 
are in antiphase with each other. This can also 
be seen in Fig. 2, which was taken with the knife 
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edge normal to the stream. Also visible here is the 
tendency for the sound waves to exhibit a pecu- 
liar “double” structure; reflections off the struc- 
ture supporting the arm carrying the edge can be 
seen crossing the lower part of the field. 


Fig. 2. Edge tone at high speed with Schlieren knife-edge 
perpendicular to jet. 

In view of the fundamental nature of this phase 
difference, it was necessary to verify that this is 
also the case in the normal low speed edge tone, 
where the density gradients are insufficiently 
large for the optical method used. To do this a 
pair of similar crystal microphones were used in 
conjunction with a pair of amplifiers and a double- 
beam cathode ray oscilloscope. By positioning 
both microphones first on one side of the jet when 
a low speed edge tone was present (e.g. at an 
efflux velocity of 9 m/s), and then on opposite 
sides, the 180° relative phase shift could be 
clearly shown. In the plane of the jet the micro- 
phone response falls to a low minimum. This mi- 
nimum was well marked, the intensity being 
otherwise more nearly constant, as in the high 
speed case (see Fig. 1 and 2), than sinusoidal. 

In order to verify that acoustic and not hydro- 
dynamic fluctuations were being measured by the 
microphone, the variation of intensity with dis- 
tance from the edge was investigated. It was 
found that the signal amplitude varied inversely 
with this distance, in accordance with the inverse 
square law for radiated energy. In addition the 
microphones used were pressure receivers so that 
if the sources are “‘simple” the pressures de- 
tected can only be acoustic, and it seems unlikely 
that sources of higher order would produce the 
simple directional pattern observed. The latter 
remark also applies to the effects of the stream 
disturbances, per se, whose pressure radiations 
are unlikely to be of comparable intensity with 
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the edge tone (even if only by virtue of the ap- 
parently higher order of acoustic source of the 
stream disturbances, i.e. a quadrupole field, 
LicuTuiy [12]). 

Thus the edge tone appears to have the struc- 
ture of a pair of equal and opposite sources, each 
radiating nearly uniformly on its side. This is so 
in both low and high speed cases; it should be 
noted that in the latter the edge is in a subsonic 
region of the stream, and that it was found that 
actual choking of the jet was incidental to changes 
in the total head of the jet producing edge tones. 
Thus so far as the actual production of the sound 
is concerned, there appears to be no fundamental 
difference between the two cases. 

That the sound should be of opposite phase on 
each side of the jet is not hard to see when it is 
realised that the motion and the sound produced 
are of the same frequency, for since the hydro- 
dynamic motions of the fluid about the wedge at 
this frequency are alternating in character, then 
so must the acoustic motions. 

But it may be asked why the equal and op- 
posite simple sources do not form a dipole with 
its characteristic figure of eight distribution, par- 
ticularly in the low speed case. 

An explanation is suggested along the follow- 
ing lines. Careful examination of the photographs 
reveals that the sound waves effectively emanate 
from a point on or near the surface of the aerofoil, 
some quarter of its chord from its leading edge. 
In a study of the low speed case, RICHARDSON [5] 
found that the maximum velocity fluctuations 
occurred adjacent to the surfaces at an appreciable 
distance from the actual edge: it may well be that 
this indicates the origin of the sound, similarly to 
the high speed case. 

Under these circumstances the surfaces of the 
wedge or aerofoil will support the fluctuating 
pressures of the flow and prevent the direct 
oscillatory movements associated with dipoles. 
Away from the edge, lateral motions* will take 
place across the plane of symmetry tending to 
cause degeneration towards dipoles. Note that 
even if such degeneration were complete, lateral 
motions would still occur. At low speeds the jet 
stream will conform to these lateral oscillations, 
but may offer some impedance at higher speeds. 

With this information, the mechanism of edge 
tones immediately suggests itself along the fol- 
lowing lines: the sound produced by the inter- 
action of the stream disturbance and the edge 


* Note added in proof. Since such motions will decrease 
more rapidly with distance than purely acoustic motions, 
they could conceivably be described as ‘hydrodynamic’ in 
nature. See CuRLE, Proc. roy. Soc. 216 [1953], 412. 
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sets up an alternating pressure gradient across 
the jet, as a consequence of the observed phase 
opposition, and in a manner analogous to that of 
“sensitive” jets produces periodic disturbances 
at the orifice which will be amplified as they pass 
downstream until they themselves produce sound 
on reaching the edge, so perpetuating the process. 
The fact that the jets of edge tones are “sensitive” 
i.e. the sound-induced disturbances undergo a 
high amplification, is of course what has been 
found in experiment, and has already been made 
use of in the associated case of the supersonic jet 
without an edge (PowELL [11}). 


3. Theoretical considerations 


Now let us consider the suggested mechanism 
in greater detail, and investigate the possibilities 
of placing the theory on a quantitative basis. 
Since the exact point from which the sound can 
be considered to emanate is not known, the arti- 
fice of an effective source distance h’ from the slit 
orifice will be introduced. It appears a priori that 
one is not entitled to assume that a certain phase 
of sound wave will initiate a stream disturbance 
of the same phase as that which gave rise to the 
sound wave in the first place; so in general let it 
be assumed that the phase difference between the 
two disturbances associated in this way with a 
sound wave is 9, the new disturbance at the slit 
naturally lagging. It will be noticed that 6 may 
be a non-integer, and greater than unity. 

A disturbance of arbitrary but identified phase 
leaving the slit will arrive at the effective source 
at atime |(1/u) dh’ later, the disturbance velocity 
u being a function of h’, and the associated phase 
of sound wave will arrive back at the slit after a 
further time h’/c has elapsed (c being the velocity 
of sound). By the definition of 9, that many 
cycles of the stream disturbance must have been 
initiated in the meantime, so that 


the frequency being f. This is a fundamental 
relationship of the theory. Since both 6 and f have 
to be found, the relationship is insufficient in it- 
self. However, it is clear that a further condition 
exists if an edge tone is to be self-maintained. 
This is that after following a complete cycle of 
the motion,—extraction of energy from the 
stream by the process of instability, its conversion 
to acoustic energy and the reaction upon the jet 
to provide new disturbances—the magnitude of 
the quantity considered shall not have changed. 
The motion may be represented functionally as 
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Fig. 3. Functional diagram of edge tone mechanism. 


shown in Fig. 3, the main components of which 
are self-explanatory. 

The hydrodynamic back reaction from the 
motion in the neighbourhood of the edge has been 
included since the happenings there must influ- 
ence the later stages of the flow, and hydrodyna- 
mic effects at the orifice may be appreciable under 
certain circumstances, when it may not be possible 
to distinguish clearly between hydrodynamic and 
acoustic disturbances in this region, e. g. when h’ is 
small and the acoustic wavelength large. (In fact 
Brown has considered this back reaction on the 
stream to be the dominant effect.) Also there will 
be some mutual interference of the sound waves 
across the stream, so reducing the effects at the 
orifice although it seems that the presence of the 
stream mitigates this effect. Thus if q,, denotes the 
amplification of stream disturbances over the dis- 
tance h’ and the 7's denote an “ efficiency”’ of the 
other processes, the criterion for the establish- 
ment and maintenance of an edge tone can be 
written as i 

qu Ns Hr Ha 21 (2) 
where the hydrodynamic reaction and the mutual 
cancellation have been omitted for simplicity: it 
is assumed that these effects are small in the 
region we shall consider, but this may not always 
be the case. Alternatively, the latter might be in- 
cluded in 7. 

The problem therefore is to obtain a solution 
to equations (1) and (2), and we shall now con- 
sider the difficulties which these present and how 
approximate solutions might be derived. In turn 
we shall consider the initiation of the stream dis- 
turbance, its career and finally the sound pro- 
duction. 


a) Initiation of the stream disturbances 

The idea that the stream disturbances of “‘sen- 
sitive’ jets are induced by aerial velocity fluc- 
tuations has come about principally as the result 
of two experimental facts: i.e. maximum sensiti- 
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vity is exhibited at velocity antinodes of standing 
waves where direct pressure effects are small, and 
also that the plane of the stream disturbance and 
vortex formation show a marked directionality 
with regard to the impinging sound waves, the 
transverse undulations of the jet lying in the 
plane of the particle vibrations of the sound 
waves. It is generally agreed by experimenters 
that the stream disturbances of sensitive jets 
originate at the nozzle. ANDRADE [13] considers 
that this happens because the fluid which has 
emerged from the orifice is more or less free to 
vibrate with the surrounding acoustic medium, 
whereas that portion of the stream yet to leave 
the orifice is constrained from transverse motions. 
Hence distortions of the stream are caused, tem- 
porally, as it leaves the orifice, and these are 
carried downstream by the flow, becoming am- 
plified if unstable. It is useful to note in the latter 
respect that “sensitivity” is only observed if 
sufficient amplification takes place, and absence 
of “sensitivity” does not necessarily mean that 
no stream disturbances are created. 

It is usual in stability calculations to consider 
the initial disturbance as a displacement or an 
imposed velocity: but in fact in this case it must 
be the pressure gradients across the jet which are 
fundamentally responsible, and so long as the 
initial disturbances are relatively small, we might 
suppose that the effects will be roughly propor- 
tional to the pressure gradient across the jet. 

The sound waves acting on sensitive jets are 
usually progressive, but even if they are of the 
standing type, the pressure gradients across the 
jet will be roughly inversely proportional to the 
wavelength. But in the case of edge tones we 
have seen that the sound waves are of opposite 
phase on each side of the stream, so instead of the 
pressure gradient being zero, as in the case of a 
progressive wave proceeding from the down- 
stream direction, it is inversely proportional to 
some dimension not a great many times in excess 
of the width of the jet. Thus edge sounds may 
well be some hundreds or thousands of times more 
effective than ordinary progressive waves of like 
intensity. Since jets may be of about the same 
sensitivity, in the ordinary sense, as the ear, it is 
not surprising that edge tone motions may exist 
when the sound produced is inaudible or that in 
other words greater intensities of progressive 
waves are necessary to produce like motions in 
an edgeless jet, as observed by Brown [4]. 

The critical nature of RicHARDSON’s experi- 
ment [5] with an edgeless jet in the proximity of 
another jet performing normal edge tone motions 
is also removed: for if the slits are in the same 
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plane, unless very close together, the mutual 
cancellation of the sound waves over the inter- 
vening distance (which could be easily demons- 
trated by microphone traverses) will greatly 
reduce their effectiveness, whereas if they are not 
in this plane there will be no tendency to cause 
disturbances in this way, but only in the ordinary 
manner of sensitive jets which as we have seen is 
far less effective, especially if directional effects 
are also present. 


b) Career of the stream disturbances - 
(i) When the velocity profile is square 


The stability of a jet, having as boundaries 
two plane velocity discontinuities, was first dis- 
cussed by Rayueicu [6]. By assuming initial 
deformations of sinusoidal form, and of small 
amplitude A, then the solution of his equation (54) 
gives the amplitude at time ¢ later for a “sinuous” 
disturbance, as 


Vt 


a = Acosh wk Vt-cosk (7 hET—*) (3) 


where x is taken along the stream 


y = tanh kl/(1 + tanh kl)?, 
k = 2/A with A the disturbance wavelength, 
V = the (uniform) stream velocity, 
1 = semi-width of the jet, which we shall alter- 
natively call 6/2. 


On the other hand if the disturbances are looked 
upon as velocity-induced a factor 2zf/ukV is 
introduced. (But we are more interested in ampli- 
fication with respect to distance rather than with 
respect to time at the moment.) The amplification 
q over a distance equal to the slit width is in 
either case 


q = exp [wkd (1+ coth kl)]} 
= exp [27S u (1+ coth kl)?] 


(4a) 
(4b) 


provided that q is not too small. Here we have 
introduced the Strouhal number 


S = fd|V. (5) 


RAYLEIGH noticed that his result was contrary 
to the observations on sensitive flames which 
possess both upper and lower frequency limits, 
whereas the degree of instability found increases 
indefinitely with frequency. Arguing that the 
instability is likely to be influenced when the dis- 
turbance wavelength and the boundary layers 
are of comparable size, RayLeicu [14] included 
viscosity effects in an arbitrary manner by as- 
suming various velocity profiles of the jet. For 
example a jet with a maximum velocity at the 
centre falling uniformly to zero on each side dis- 
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plays a maximum instability for sinuous dis- 
turbances when their wavelength is about 2.56, 
decreasing with shortening wavelength until at 
A= 1.786 actual stability is found. Although this 
is a hypothetical result it does give a hint on the 
behaviour of a real jet after development of the 
boundaries due to viscosity has taken place. 

However, the first result obtained by RAYLEIGH 
quoted above can be used over that part of the jet 
close to the exit, ANDRADE’s work [15], [16] having 
shown that in the majority of sensitive jet and 
edge tone experiments the velocity profile might 
be considered as being square. On taking, for 
example a frequency of 126 c/s, and a slit width 
of 1mm, stream velocities of 380,150 and 100 cm/s 
result in disturbance amplifications (over a dis- 
tance 506) of 30.8, 8400, and 54000 respectively. 
This increased instability for the slower stream, 
at first sight perhaps surprising, can be clearly 
seen if reference is made to Plate 2 of Brown’s 
paper [17] where the vortex formation approaches 
the orifice as the speed is decreased. It is interest- 
ing to note that such rates of amplification are 
capable of developing acoustic disturbances to a 
visible size over a distance comparable to that 
observed experimentally. 


(ii) Developed velocity profile 


Some distance from the orifice the boundary 
layers will have developed to an appreciable ex- 
tent, finally merging with each other to form ap- 
proximately a steady laminar jet like that which 
issues from a line source, as calculated by 
Bick.ey [18]. By applying the method of small 
disturbances to such a velocity profile, taken as 
one-dimensional, Savic [19] found the nature of 
the neutrally stable disturbances of such a flow. 
These took the form of an alternate arrangement 
of areas of like vorticity superficially resembling 
a KArMAN vortex street, with a spacing 


A=a)b (6) 


and moving with the velocity of the point of in- 
flection of the basic velocity profile, 
uw. 2 


Sk ae 1, 
7 = ze = 1016S". (7) 


The “‘constants”’ a and 6 are obtained from the 
basic flow, being 


/ K2\ "l/s Ms 
a = 0.4543 ra > b6=0.2751 % (8) 
yx yx 


where K is the momentum of the jet per unit slit 
length per unit time in the x direction, and y the 
kinematic viscosity. As Savic pointed out, if the 
exit velocity profile can be considered square, 
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then K = V*6 where V is the exit velocity, hence 
equation (7). 
Despite the fact that this result pertained to 


neutral disturbances, assumed small, Savic pro-| 


ceeded to compare the disturbance wavelength so 
found with those observed by Brown; and he 
obtained good agreement. It would seem, there- 
fore, that the neutral solution is not far different 
from the highly unstable motions of sensitive jets 
and edge tones, so far as disturbance velocity is 
concerned. 

Now let us examine the implications of the 
solution obtained by Savic, who considered that 
the stability of the jet broke down at a certain 
point depending on the frequency. Except for 
very high frequencies the flow is always unstable 
on leaving the orifice, and the disturbances, if 
still “‘small”’, will sooner or later enter a stable 
regime; if they are large stability may never be 
achieved although the degree of instability need 
not necessarily be great. That the boundary found 
by Savic is a transition to stability can be seen 
by reflecting that since the stream is decelerating 
the wavelengths corresponding to a given fre- 
quency may decrease, while at the same time the 
layers of vorticity are thickening, so that the 
trend must be towards stability. Savic’s solution 
strictly refers to a single section of the jet, namely 


x/d = 0.02650-R/S (9) 
where R is the Reynolds number V6/y, x being 


measured from the hypothetical line source. 
ANDRADE [16] showed beyond a certain distance 
BIcKLEY’s solution could be applied to a real jet 
provided the hypothetical source was taken a dis- 
tance 0.65 Rd behind the real slit. Thus if this 
stability boundary is not to be “‘behind”’ the real 
slit 


S < 0.0407. (10) 


Brown worked in the ranges 0.0015 <S <0.5 for 
sensitive jets and 0.015 <S <0.13 for edge tones. 
It is interesting to see therefore that had the dis- 
turbances remained small (and even if ANDRADE’s 
estimate of 0.65 Rd is not very accurate), in 
many cases they would have become stable on 
reaching the region of steady laminar flow, if not 
during the transition. Brown [17] in fact reported 
the subsidence of symmetrical disturbances, 
where the trend to stability is greater. The im- 
portance of both Reynolds number and Strouhal 
number is now apparent, and failure to recognize 
this has resulted in various apparent disparities 
in the past. Since the nature of the stability will 
also be dependent upon the magnitude of the 
disturbances, a further parameter linking the 
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disturbing sound with a stream parameter should 
also be included for strict dynamical similarity. 
However, it is instructive to make a comparison 
on an R-S basis, when, taking BRown’s measure- 
ments of relative vortex velocity for example 
(Fig. 8 Brown), a substantial degree of dynami- 
cal similarity will be found over a frequency 
increase of fourfold. 

We can now explain the limiting frequencies for 
“sensitivity”’: for if the frequency is too high 
stability will set in at an early stage, whereas if 
it is too low the rate of amplification will be too 
small. These frequency limits, with a region of 
maximum overall instability between them, are 
of great importance with regard to the solution of 
equation (2). Few data are readily available from 
which to find the rate of amplification and dis- 
turbance velocity for the changing velocity pro- 
file, but equation (7) gives a reasonable approxi- 
mation under certain circumstances. 


c) The sound production 


A quantitative estimate of the sound produc- 
tion is required to enable a complete theoretical 
solution to the basic equations to be made. Al- 
though a considerable advance has been made by 
LiGHTHILL in the estimation of aerodynamic 
sound, such an estimate is not yet possible by 
virtue of lack of detailed knowledge of the fluid 
motions about the edge, since “large” disturb- 
ances and the influence of the edge itself must be 
taken into account quantitatively. An analysis of 
these would locate the effective source and par- 
tially determine the phase difference 0, which we 
shall for convenience write 


6=N+p (11) 


where NV is an integer and p< 1. The value of p 
will be essentially constant, although it may vary 
somewhat with conditions at the edge, since the 
stream here may be little more than sinuous, with 
vortex development delayed until adjacent to the 
faces of the wedge or at the other extreme early 
vortex growth may be commencing to disinte- 
grate. The assumption that a maximum sound 
pressure is generated just as a vortex reaches the 
edge leads to p=1/, since this would lead to a 
maximum stream deflection to the opposite side. 
The value of p so found is in satisfactory agree- 
ment with Brown’s observations, the non-integral 
part of h/A always lying between 0 and 0.5 except 
for one value (0.9). Also, lacking detailed infor- 
mation we might assume that the effective source 
is at the edge itself, distant A from the slit, i.e. 


h’=h. (12) 
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d) On the form of the solution 


It is clear that further work remains to be done 
before the general solution can be obtained, but 
nevertheless several features peculiar to edge 
tones may be deduced. Since the amplification q 
is more frequency-dependent than the other fac- 
tors of equation (2), it is clear that it will largely 
determine the parameters of the motion which 
must take place near the frequency for maxi- 
mum amplification. This region roughly coincides 
with that over which sensitivity is found. But the 
phase relationship (1) must be satisfied, and since 
p is practically constant, it follows that N may 
have one or more values. N clearly indicates the 
number of complete disturbance wavelength pre- 
sent between the slit and the edge in any stage, 
and it so happens that it is numerically equal to 
the stage number according to BRown’s conven- 
tion [4] as may be seeen from his photographs. 


The possibility of more than one edge tone 
being produced at a time can be realised, since N 
may have more than one acceptable value. Also 
some insight into the mechanism of the “jumps” 
in frequency may be gained. If one of the para- 
meters be varied, q will move away from the 
optimum value for equation (2) since 6 will be 
fixed, but gradually conditions will change so 
that the neighbouring value of § will become pre- 
ferable. But the basic flow for this new mode is 
the steady jet flow together with the existing 
mode. If the latter is “large” it will introduce 
additional damping of the new mode so tending 
to suppress it. (This is easily visualised if the 
original mode develops into the vortex state.) But 
this state of affairs cannot persist, since the am- 
plitude of the original mode will decrease as 9 and 
q move away from the optimum values, so that 
at some point the weaker mode will meet the gain 
conditions of (2). The original mode will then be- 
come damped as the new mode establishes itself 
and will now dominate, i.e. a jump has taken 
place. Clearly the jump in the reverse direction 
will be similarly delayed, so producing the cha- 
racteristic hysteresis. It is easy to see, when 
nearing the end of a hysteresis region, how a 
sharp noise may suddenly precipitate a change in 
note. Thus the immediate past history of the 
motion has some influence on N. 

Since the value of N can be obtained from the 
condition of nearness to maximum “sensitivity” 
—or directly from Brown’s observed stages— 
and the value of p might be taken as !/,, we can 
obtain an approximate solution, applicable under 
restricted conditions although a general theore- 
tical solution cannot yet be made. 
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Let us assume that h=h’. We cannot therefore 
expect good agreement with experiment if h is 
too small since then the approximation will be 
poor, especially if the vortex development is 
delayed. 

If the Reynolds number is moderately high, 
and h not too large, we might take RAYLEIGH’s 
solution for disturbance velocity to give a value 
not far different from the mean value. Taking 
this from equation (3) and placing it in equation 
(1) together with the other assumed values gives 


WeVOEN P25 mee 
Sap aa 1+ cothkl 


(13) 


where u is that mean velocity, and with the ap- 
proximate value from (1) 


kl = x (N + 0.25) d/h. 


Similarly taking Savic’s solution (7) as giving 
approximately the mean velocity, i.e. with the 
Reynolds number and 6/h not being too large 


V [N+0.25]* /6\% 

LEM A amare (i) 

A comparison can be made with Brown’s 

experimental results, these being shown in Tables 

I and IT. In view of the restrictions, equation (13) 
has been limited to the region 6/h=0.1. 


It is seen, first of all, that the estimated fre- 
quencies are of the right order, the errors being 


(14) 


Table I 
d6=0.1 cm 


f estimated 
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Table II 
6=0.1 cm, h=3.5 em, eq. (21) 


estimated 


observed 


considered reasonable in view of the approxima- 
tions made. Also the magnitude of the jumps are 
in good agreement. In the region 6/h=0.1, the 
formula (13) gives the better result, as might be 
_ expected from the conditions attached to the 
assumptions. 

It is interesting to note that agreement with 
RIcHARDSON’s results [5] is poor: for here h was 
small, and the maximum velocity fluctuations 
were found to be an appreciable distance from 
the edge, the effective acoustic source being pre- 
sumably similarly placed. 

These formulae cannot be applied to the high 
speed case as photographed. This is because the 
flow is turbulent, and insufficient data have so far 
been collected to enable a determination of the 
mean disturbance velocity to be found. It is 
interesting to note that as the velocity is raised 
from an ordinary low-speed edge tone value, after 
several jumps in frequency the note quite quickly 
becomes rough so that no discrete note may be 
distinguished as general turbulence ensues. But 
if the jet velocity is still further increased a new 
edge tone appears, of the type photographed. 
(The jet need not have reached sonic speed.) 
Owing to the additional damping due to the 
“eddy” viscosity, amplifications are compara- 
tively small, the sound production being more 
efficient by virtue of the high stream velocities 
and higher frequencies so that the gain criterion 
(2) may be met. The very high intensity of the 
sound confirms this idea. The mechanism at the 
orifice undergoes a change if sonic efflux ve- 


A.POWELL: ON EDGE TONES AND ASSOCIATED PHENOMENA 


241 


locity is attained. Instead of the jet moving 
bodily relative to the nozzle, crudely speaking, as 
in the low speed case, the initial disturbances are 
now generated by the acoustic pressure fluctu- 
ations at the nozzle causing the angles at which 
the jet boundaries leave the orifice to fluctuate. 


4, Associated phenomena 


Under this heading we shall discuss very briefly 
some phenomena which appear to be related to 
edge tones in that conditions analogous to equa- 
tions (1) and (2) are applicable. 


(i) The case of the edgeless supersonic jet 
(PowE LL [11]) mentioned in paragraph 2. Here 
the effective transmission of the sound to the 
orifice is dependent upon the interference be- 
tween the individual sources located near the 
regularly spaced shocks. Hence here 7; is acutely 
frequency dependent (for given other conditions) 
and the motion might be expected to take place 
near a maximum of 7. This is what was found, 
N taking the values 2, 3 in the cases examined 
with p effectively about 7/,. 


(ii) “‘ Jet tones”. Here no edge is present, and 
the position of the effective source, which is de- 
pendent upon the magnitude of the disturbances, 
is probably located near the point at which vortex 
development proper commences. Since the effi- 
ciency of the sound production is low, the overall 
amplitude of the motion is comparatively small, 
resulting in a weak radiated sound. Under these 
conditions the back reaction or “‘ wake disturb- 
ance”’ (Fig. 3), as envisaged by Brown [2] may 
be of greater importance. Also in this class is the 
whistling of jets when a breakaway at the exit is 
present, this being susceptible to alternating pres- 
sure gradients resulting from the flow. The 
same may be true of the periodic wakes formed 
behind bluff bodies. It is interesting to note that 
these motions may all be eliminated, as well as 
those of edge tones, by introducing additional 
turbulence into the stream, or causing premature 
turbulent mixing so that qn is reduced below the 
critical value demanded by equation (2). 


(iii) The influence of local resonators on an 
edge tone is to artificially modify 74, the trans- 
mission of sound to the orifice. It is perhaps 
better to write 


gh *s*(He+ Hr) ya 21 
where 7 is the additional part due to reflectors 
or resonators. If 7, is small and positive, one 
would expect a larger amplitude of the motion 
and slightly increased sound emission, and an 
increased tendency to jump to the next mode. 
This is easily observed with a low-frequency edge 
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tone by placing a reflector close to the jet, al- 
though far enough away not to interfere hydro- 
dynamically. On the other hand when 7, is an 
alternating function, easily obtained for the high 
speed case, and @ has a limited possible range (as 
appears to be the case) then for certain reflector 
positions the gain criterion just stated can be met, 
but not so over the intervening distances. This 
may be shown very strikingly by moving the 
hand to and for as a reflector, there being a series 
of regions for which no sound is generated and 
between these the sound builds up to a maximum 
and then falls away again as the intermediate spa- 
ces are traversed. The writer considers that the 
effects just described are convincing evidence of 
the importance of acoustical effects in the edge 
tone mechanism, at least in the high speed case 
and over the investigated range for the low speed 
case. 

The final case we shall consider is that of a 
highly tuned local resonator, that is in which 7, 
has at least one high maxima. Here the gain cri- 
terion may be satisfied, by virtue of the size of 7;, 
even though q may be small and far removed 
from the region of “visible” sensitivity. In the 
case investigated by the writer the placing of 
reflectors on each side of the jet resulted in a 
single note appearing over a range of low speeds, 
fading out to reappear at a much higher jet speed 
not necessarily above the choking speed. The 
latter case is shown in Fig. 4, where the large dis- 
turbance amplitude developed will be noted, the 
jet debouching at sonic speed. 

A similar state of affairs seems to exist in an 
organ pipe, where the resonant frequency may be 
much less than the corresponding edge tone value 
(Mercer [20]); the effect of the resonator then 


Fig. 4. Parallel reflectors resulting in enhanced motion at 
fixed frequency. 
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greatly exceeds that of the sound proceeding 
directly from the motion at the upper lip of the 
pipe, so forcing the frequency. However, during 
the starting process, relatively strong harmonics 
may be developed which may later diminish. 
These may arise because their amplification in the 
stream is greater (being nearer the region of nor- 
mal sensitivity) so tending to set up these modes, 
but the building up of the resonance at the funda- 
mental frequency may partially damp these out 
rather analogously to the jumps of edge tones. 
The “nicking” of the orifice to remove unwanted 
higher frequency noise, possibly arising from mis- 
cellaneous resonances may well result in the 
damping of the associated disturbances (by “eddy 
viscosity’) while hardly affecting the lower fre- 
quency disturbances of the fundamental. 

But we shall not trespass further into these 
interesting fields: they are studies in themselves. 
It is hoped that the foregoing may perhaps pro- 
vide a new facet to these problems on which 
much work remains to be done. 


5. Concluding remarks 


The observation that the sound produced by 
edge tones is of opposite phase on the two sides 
of the jet makes possible an explanation of the 
phenomenon in terms of the sensitivity to sound 
of the jet, it being particularly susceptible to this 
form of acoustic disturbance. This is because of 
the relatively large pressure gradients existing 
across the jet, the mechanism being in contrast 
to existing current theories which are essentially 
hydrodynamic. Such hydrodynamic effects per- 
haps causing “sensitivity” directly or reacting 
upon the stream may be important under certain 
conditions, for example when the edge is very 
close to the exit. Although a method by which a 
general solution may be effected has been indi- 
cated, it is not applicable in practice until each 
part of the mechanism can be placed on a quanti- 
tative basis, the motion of the fluid in the neigh- 
bourhood of the edge and the associated sound 
production in particular still presenting a major 
problem. However even in its present form the 
theory is capable of explaining the salient charae- 
teristics of edge tones and the agreement with 
experiment of the approximate solutions are 
sufficiently encouraging to warrant further in- 
vestigation, perhaps along the lines indicated, or 
of the nature of a “crucial test”. With regard to 
the latter, the previous major criticisms of this 
type of mechanism have been explained, the 
reason for the negative result of RicHAaRDSON’s 
experiment also being demonstrated experimen- 
tally. It is more difficult to give positive experi- 


ACUSTICA 
Vol. 3 (1953) 


mental proof that an acoustic linkage is normally 
involved, although the susceptibility to the pre- 
sence of reflectors not interfering with the stream 
hydrodynamically is direct evidence of the im- 
portance of acoustical effects. 
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SCHALLABSTRAHLUNG VON PLATTEN, 
DIE ZU BIEGESCHWINGUNGEN ANGEREGT SIND 
Von K. GOSELE 


Institut fiir Technische Physik, Stuttgart 


Zusammenfassung 


Die Schallabstrahlung einer zu freien Biegewellen erregten Platte wird, wie bekannt, von dem 
Verhialtnis der Biegewellenlange Jz zur Wellenlange A; des Luftschalls gleicher Frequenz beein- 
fluBt. Es wird rechnerisch die Abhiangigkeit der Schallabstrahlung einer biegewellenerregten 
Platte von deren Abmessungen behandelt. Dabei wird gezeigt, das kleine Abmessungen zu einer 
relativ starken Abstrahlung auch fiir die Biegewellen fiihren, deren Wellenlainge kleiner als /z ist. 


Summary 


The radiation of a plate excited in free flexural vibration is influenced by the ratio of the 
wave-length / of these waves to the wave-length Ay of the sound in air. The relation between 
the radiation of a plate in flexural vibration and its dimensions is calculated, and it is shown 
that small dimensions lead to a relatively strong radiation even for flexural waves whose wave- 
length is less than A. 


Sommaire 


On sait que l’émission sonore d’une plaque excitée suivant des oscillations libres de flexion 
dépend du rapport de la longueur des ondes de flexion Ag 4 la longueur d’onde Ay, du son aérien 
de méme fréquence. On calcule dans le présent article l’influence des dimensions de la plaque sur 
l’émission sonore d’une plaque oscillant a la flexion. On démontre que |’émission est rela- 
tivement intense lorsque la plaque a de petites dimensions, méme dans le cas d’ondes de flexion 


dont la longueur d’onde est plus petite que /,. 


1. Einleitung 


Die Schallabstrahlung von Platten, die Biege- 
schwingungen ausfiihren, ist bisher verhaltnis- 
mabig wenig behandelt worden. 

Beobachtungen der letzten Jahre itiber die 
Schalldammung von zweischaligen Decken- und 
Wandkonstruktionen [5], [7], [8] haben jedoch 
gezeigt, daf} diese Frage eine erhebliche prak- 
tische Bedeutung hat. Vom _ bautechnischen 
Standpunkt ist es haufig erwiinscht, daB zwischen 


den Schalen (1) und (2) einer zweischaligen Decke 
oder Wand (siehe Abb. 1) an einzelnen Punkten 
eine feste Verbindung (3) besteht. 


(2) 
Abb. 1. Aufbau einer Decke mit zweiter Schale. 


24.4. 


Eine solche Verbindung ist nach neueren Un- 
tersuchungen vom akustischen Standpunkt [5], 
[8] dann zulassig, wenn die eine der beiden Scha- 
len eine geringe Biegesteifigkeit aufweist. Die von 
der Schale (1) auf die Schale (2) tibertragenen 
Biegeschwingungen fiihren in diesem Fall zu einer 
relativ geringen Luftschallabstrahlung. Die un- 
geniigende Schalld4mmung von Decken und 
Wanden kann durch eine derartige Verkleidung 
mit abstrahlhemmenden, biegeweichen Schalen 
wesentlich verbessert werden. 

Die Abstrahleigenschaften von Platten, die 
Biegeschwingungen ausfiihren, sind qualitativ 
bekannt. Sobald die Wellenlange Ap der Biege- 
wellen kleiner als die Wellenlange Ax des Luft- 
schalls gleicher Frequenz ist, wird die Abstrah- 
lung wesentlich geringer als normal (siehe [1], 
[4]). Quantitative Unterlagen iiber das Abstrahl- 
verhalten fehlen jedoch noch. Im folgenden soll 
ein Beitrag zur Berechnung der abgestrahlten 
Luftschallenergie einer zu freien Biegewellen er- 
regten Platte gegeben werden. 


2. Unendlich ausgedehnte Platte 


Die Schallabstrahlung einer unendlich ausge- 
dehnten ebenen Platte, die Biegeschwingungen 
(Schnelle- Amplitude vy, Biegewellenlange /z) aus- 
fiihrt, ist von L. CREMER [6] und spater von 
J. Brittovuin [8] berechnet worden. Die abge- 
strahlte Schall-Leistung N der Flache S der Platte 


ergibt sich danach zu: 


Nae vs Esust tog fir Ap 2dr, (1) 
2 V1 —(Az/Az)? 
N=0 fir Ap<dzt. (2) 


Zur anschaulichen Darstellung soll im folgen- 
den ein dimensionsloser Abstrahlfaktor s einge- 
fiihrt werden, der die abgestrahlte Leistung N, 
die von dem Flachenstiick S einer unendlich 
groBen ebenen Platte abgestrahlt wird, in Be- 
ziehung setzt zu der Leistung No, die von S ab- 
gestrahlt werden wiirde, wenn alle Plattenpunkte 
konphas schwingen (siehe Abb. 2). Im letzten 
Falle strahlt die Platte eine ebene Welle ab, deren 
Fortpflanzung in Richtung der Plattennormalen 
erfolgt. Die Schwingungsschnelle v, sei bei beiden 
Platten dieselbe. 


Dann gilt: 
No = . vo: S, (3) 
s =1///1—(Ax/As)? fir As = Az, (4) 
s=0 fiir Ap < AL. (5) 
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Abb. 3 zeigt den Verlauf von s in Abhangigkeit 
von Ap/Ar. Die Abstrahlung hért fiir Schwin- 
gungen, deren Wellenlange Az < Az ist, vollig auf. 
Fiir Az = Ar erreicht sie einen Maximalwert, um 
mit héher werdendem Verhaltnis Az/A, wieder 
abzusinken. Die bisherigen experimentellen Be- 
obachtungen an endlichen Platten zeigen nun, 
daB auch fiir Ax < Az Luftschall abgestrahlt wird. 


Abb. 2. Zur Definition des Abstrahlfaktors s = N/ Np. 


Abb. 3. Verlauf der Schallabstrahlung bzw. des Abstrahl- 
faktors s bei einer unendlich groBen, zu Biege- 
schwingungen (Ag) erregten Platte. 


Als Ursache fiir die Abweichung zwischen dieser 
Beobachtung und dem idealisierten Fall der 
Theorie kommen vor allem die endlichen Ab- 
messungen der Platten in Frage. Die obige Defi- 
nition des Abstrahlfaktors soll fiir diesen Fall et- 
was erweitert werden. N soll dabei die abge- 
strahlte Schalleistung der gesamten, endlich 
groBen Platte von der Gesamtflache S sein. Sie 
soll wieder mit der oben angegebenen Leistung N, 
eines gleich groBen Flachenstiickes einer unend- 
lich groBen, konphas schwingenden Platte ver- 
glichen werden: 


s = N/N,. (6) 
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3. Platten mit endlichen Abmessungen 


a) Fortschreitende Wellen 


Zur Berechnung des Abstrahlverhaltens einer 
Platte mit endlichen Abmessungen soll eine 
rechteckige Platte (Kantenlangen a, b) angenom- 
men werden, iiber die eine fortschreitende Biege- 
welle in Richtung y wandern mége. Fortschrei- 
tende statt der im praktischen Fall auftretenden 
stehenden Wellen wurden zunachst gewahlt, um 
von den verschieden wahlbaren Randbedingun- 
gen an den beiden Plattenenden frei zu kommen, 
die das Ergebnis beeinflussen. Die benutzte Rech- 
nungsannahme kann man sich durch eine unend- 
lich groBe Platte realisiert denken, die zu fort- 
schreitenden Biegewellen angeregt ist. Diese 
Platte werde durch eine zweite, starre Platte ab- 
gedeckt, die mit einer rechteckigen Offnung mit 
den Kantenlangen a und 6b versehen sei, so dab 
nur iiber diese Offnung eine Abstrahlung méglich 
ist. 

Die Biegeschwingungen der Platte P mégen die 
Schnelle v, aufweisen: 


i(wt—kpy) 
ede (7) 


Up = Vy" e 


x 
Abb. 4. Zur Bestimmung der Schallabstrahlung einer end- 
lich groBen Platte mit Biegeschwingungen v p. 


Die von der Platte in den Halbraum abge- 
strahlte Schalleistung N soll bestimmt werden. 
Zur Lésung wird, wie es Abb. 4 zeigt, fiir einen 
im Halbraum befindlichen beliebigen Aufpunkt O 
mit den Koordinaten R, x und f das Geschwin- 
digkeitspotential ® des Schallfeldes bestimmt, 
das von der Platte herriihrt. Die Berechnung von 
® ergibt (siehe z. B. [2]): 


i(wt—k, R) 
-e Sah 


(8) 
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wobei A und B folgende Abkiirzungen bedeuten: 


2 sin [kz (a/2) cos x] 

Ae ki cos x 3 (°) 
2 sin [(b/2) (kx cos B—ka)| 

a ki cos 8 —kp O. 


und wobei RS a bzw. 6 sein muB. 

Die von der Platte abgestrahlte Schalleistung 
errechnet sich aus der Leistung N, die durch die 
Halbkugelflache vom Radius R_ hindurchwan- 
dert. Sie ergibt sich zu: 


c= Bp=x 
1 ae 
N= gpacekins { 4°ds | Begg, (11) 
x=0 p=0 


Die Integration iiber den Winkel « 1a8t sich 
nach STENZEL [2] durchfiihren, wobei sich fir 
a>dAr der Wert za/kx ergibt. Unter Einfiihrung 
des oben definierten Abstrahlfaktors s wird (11) 
zu: 


ks * ein? [b/2) (kz. cos 8 —ks)] 
K bhp (1 —kz/kp cos B)? 


0 


dg. (12) 


Statt ky und kz sollen im folgenden die Frequenz 
f, die Grenzfrequenz f,, der Platte (nach CREMER 
[3]) und fiir b eine kennzeichnende Frequenz f, 
eingefiihrt werden: 


ke =" VF-fe (13) (14) 


wobei die Fortpflanzungsgeschwindigkeit cp der 
Biegewellen nach CREMER [3] in folgender Weise 
von der Frequenz abhangig ist! 


rena 05 
Damit wird (12) zu: gtrenwe (16) 
tf hee [aioe (2 cos B 1)| 5 


ed fm) 


In der Abb. 5 ist der Integrand von Gl. (16), 
einschlieBlich des vor dem Integral stehenden 
Gliedes, in Abhangigkeit von f dargestellt, was bis 
auf einen konstanten Faktor der Winkelverteilung 
der abgestrahlten Schalleistung entspricht. Die 
Werte sind fiir eine groBe Platte (fy/fy bzw. 
b/Agr = 30) und fiir eine kleinere Platte (fz, / fy baw. 
b/Azr = 3) fiir verschiedene Frequenzen angegeben; 

1 Es ist dabei von dem Fall ausgegangen, daf} eine Platte 


zu freien Biegeschwingungen angeregt werde, wobei die 
Abstrahlung fiir diese Schwingungen interessiere. 
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Agrbedeutet dabei die Wellenlange bei der Grenz- 
frequenz f,,. Die Grenzfrequenzen beider Platten 
sind gleich groB gewahlt. Zunachst ist der Unter- 
schied in der Schallabstrahlung offensichtlich, so- 
bald f von Werten f> fz zu solchen f < fy 
wechselt. Die stark gerichtete Strahlung bei f > fir 
verteilt sich fiir f < fz einigermafen gleichmaBig 
auf die verschiedenen Winkelbereiche. 
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1. Fir Platten, deren f<f,, ist, vereinfacht 
sich Gl. (16) fiir fy <f<fg zu: 


1 fo or 
Tacainfle u eee i 


2. Fir die Grenzfrequenz (f=f,) ergibt eine 


Naherungsrechnung: 


(17) 


s = 0,94 Vfer/ fo: (18) 


60 60 
[ Tf gr= 2 
b/ gr= 30 


40-F 40 


Zoi 20K 


m/2 1 


f/f gr= 2 


b/degn=3 3. Fir Frequenzen geniigend ober- 


halb der Grenzfrequenz wird s: 


s=1/V1—(fe/f), 


d. h. der Abstrahlfaktor nimmt ober- 
halb der Grenzfrequenz denselben 


Wert an wie fiir die unendlich groBe 


(19) 


20+ b/Rgn = 30 20+ 


m/2 T 


Platte. Dabei ist im Fall 2 und 3 vor- 
ausgesetzt, daB f) < fe ist. 

In der Abb. 6 ist der Verlauf des 
Abstrahlfaktors in Abhangigkeit von 
mt f/fer baw. Ap/Ax fir drei verschieden 


: f/f qr=05 
b/Agr= 30 


grobe Platten angegeben. Der EinfluB 
der Plattenabmessung auf die Ab- 
strahlung ist gut ersichtlich. Je klei- 
ner die Plattenabmessungen sind, um 
so gréBer ist die Abstrahlung im Be- 


f/f gr =05 
b/Agn= 3 


7/2 Mg 


eo rere f <= fer. 


b) Stehende Wellen 


Die oben skizzierte Rechnung 1aBt 
sich auch fiir stehende Wellen aus- 


m/2 54 n/2 


fiihren, wobei folgender Ansatz fiir die 
x Biegeschwingungen der Platte P ge- 
wahlt sei: 


Abb. 5. Winkelverteilung der abgestrahlten Schallenergie 


fiir eine endlich groBe, zu fortschreitenden freien l 
Biegewellen erregten Platte fiir verschiedene Fre- 
quenzen f. Als Ordinate ist der Integrand von 


B”) 


oun : (20) 


Gl. (16) einschlieBlich des dort vor dem Integral 


stehenden Gliedes aufgetragen; 


links: Platte mit grofen Abmessungen (b/A,,=30), 
rechts: Platte mit kleinen Abmessungen (b/A,,=3). 


Ein Vergleich der Werte fiir die beiden ver- 
schieden groBen Platten in Abb. 5 zeigt den Ein- 
fluB der Abmessungen der abstrahlenden Platte 
auf die abgestrahlte Energie. Die starkere Ab- 
strahlung im Gebiet f < fz bei den kleineren Ab- 
messungen ist leicht zu erkennen. (Man achte auf 
die verschiedenen Mafistabe der Ordinaten- 
achsen.) 


Naherungswerte fiir s 


Der Wert des Integrals (16) 1aBt sich fiir ein- 


zelne Teilbereiche von f naherungsweise angeben: 


10 \ogs 


01 1 10 
— for 
Abb. 6. Abhangigkeit des Abstrahlfaktors s einer Platte 
endlicher GréBe, die zu freien, fortschreitenden 
Biegewellen angeregt ist, von f/f,, bzw. Ap/A, und 
von der PlattengréBe b. 
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das Ergebnis ist: 


x{_.. (6 
Re 4 a (5 ki cos Bp + ks) 
ab (ki cos 6 + kp) 
: (21) 


sin ih ks | 
S | ala | ap: 


Dabei ist die Abstrahlung nicht auf eine konphas 
schwingende Platte von der Schnelleamplitude 


Ss 


1 é 
Umax, sondern von Ses bezogen. Es ist das 


energetische Mittel der bei der stehenden Welle 
raumlich zwischen 0 und vmax schwankenden 
Schnelleamplitude zum Vergleich herangezogen. 

Fir die weitere Rechnung, insbesondere fiir 
Naherungswerte sind im folgenden zwei Spezial- 
falle ausgewahlt worden. Die Plattenschwingung 
habe an den beiden Plattenenden 

1. einen Schwingungsknoten, 

2. einen Schwingungsbauch. 


Die zugehérigen Lésungen ergeben sich zu: 
Lésung 1 (Schwingungsknoten am Plattenrand) 


2 fo f sin®[x(fife) e058] 4 
fe J (Fifer) 005% 


B. (22) 


Naherungsweise: 


w_ 1) hk C—flfe) 
2a fer (1 — f fr) 


s=1/VI—f,/f fir f> fer > fo: 


Lésung 2 (Schwingungsbauch am Plattenrand) 


fir fy<f<fer, 


(23) 


4 


2 ff [e082 pcos? [x(f/f,) €05 5] 
“aE fe J —Cifrcot se 4) 
Naherungsweise: 
s= pF fo—fifp)"* fir fy < Jf <fon 
Jer (25) 


s=1/Vl—f,/f fiir f >for > fo. 

In Abb. 7 sind die Werte fiir eine Platte mit 
stehenden und mit fortschreitenden Biegewellen 
einander gegeniibergestellt. Der Verlauf von s ist 
fiir fortschreitende Wellen und fiir stehende 
Wellen mit Schwingungsknoten am Plattenende 
(Lésung 1) bis auf eine gewisse Verschiebung 
unterhalb der Grenzfrequenz gleich. Stehende 
Wellen mit einem Schwingungsbauch an den 
Plattenenden ergeben dagegen unterhalb der 
Grenzfrequenz eine wesentlich geringere Ab- 
strahlung als fortschreitende Wellen. 
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Fee LV fgr 
Abb. 7. Vergleich des Abstrahlverhaltens einer Platte end- 
licher GroBe (b/A,, = 30) bei stehenden und bei fort- 
schreitenden freien Biegewellen; 


Kurve (1): stehende Biegewellen (Schwingungs- 
knoten am Plattenrand), 


Kurve (2): stehende Biegewellen (Schwingungs- 
bauch am Plattenrand), 


Kurve (3): fortschreitende Biegewellen. 


4, Diskussion der Ergebnisse 


a) Verhalten fiir f> fer 


Oberhalb der Grenzfrequenz, bzw. fiir Az > Az, 
hat die GréBe einer abstrahlenden Wand- oder 
Deckenflache nahezu keinen EinfluB auf die 
GréBe der abgestrahlten Energie, solange f> fy 
ist. Der Abstrahlfaktor s ist in diesem Gebiet nur 
wenig frequenzabhangig. Fiir hdhere Frequenzen 
fallt er etwas ab. Die Annahme von L. CREMER 
[6], wonach die Abstrahlung mit wachsender 
GréBe von Ag/Ar zunehmen soll, trifft offenbar 
fir f> fg nicht zu. |] 


b) Verhalten fiir f = fer 


In der Nahe der Grenzfrequenz weist der Ab- 
strahlfaktor s ein Maximum auf, das um so héher 
ausgebildet ist, je gréBer das Verhaltnis f;,/f, 
baw. b/Agr ist. 

Je gréBer die abstrahlende Flache b ist, um so 
gréBer ist die abgestrahlte Leistung je Flachen- 
einheit. 


c) Verhalten fiir f <fer 


Unterhalb der Grenzfrequenz sinkt die Schall- 
abstrahlung stark ab. 

Sie nimmt jedoch iiberraschenderweise nicht 
mit abnehmender Frequenz stetig ab. Vielmehr 
strebt sie einem frequenzunabhangigen Grenz- 
wert zu, der fiir fortschreitende Wellen folgende 


GroéBe hat: 
(26) 


d. h. der Abstrahlfaktor hangt nur von dem Ver- 
haltnis fo/fz- oder, anders ausgedriickt, von dem 
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Verhaltnis der Wellenlange A, bei der Grenz- 
frequenz zur Abmessung b der abstrahlenden 
Flache ab. Je gréBer die Ausdehnung b der Flache 
ist, um so geringer ist der Abstrahlfaktor. Dies 
fihrt zu dem eigenartigen Ergebnis, da unter- 
halb der Grenzfrequenz eine Vergrésserung von 
b und damit eine Vergré®erung der Plattenflache 
bei gleicher Schwingungsamplitude zu keiner Er- 
héhung der in den Luftraum abgestrahlten 
Schallenenergie fiihrt. 

Die vorstehenden Uberlegungen gelten jeweils 
nur fiir den Frequenzbereich f> fy. 


d) Praktische Folgerungen 


Die Grée der in der Bauakustik auftretenden 
Abstrahlminderung liegt nach den vorstehenden 
Uberlegungen giinstigenfalls in der Gré®enord- 
nung von 20 dB. Sie wird allerdings im prakti- 
schen Fall durch zwei Effekte verringert werden: 


1. Durch meist vorhandene und nétige Verstei- 
fungen (z. B. Holzleisten) bei den in Frage kom- 
menden diinnen Schalen. 


2. Beeinflussung der freien Schwingungen einer 
abstrahlhemmenden Schale durch einzelne feste 
Verbindungen mit einer zweiten Schale. Auf 
diesen EinfluB hat J. BriLLovurn [8] hingewiesen. 


Die abstrahlhemmende Wirkung ist nur bei 
Platten von praktischer Bedeutung, die eine ge- 
nugend hohe Grenzfrequenz aufweisen. Abb. 8 
mige dies verdeutlichen. Dort ist eine Wandflache 
von 3 x 3m? angenommen, die aus verschiedenen 
Baustoffen bestehen soll. Der Abstrahlfaktor s 
ist in Abb. 8 fiir diese verschiedenen W ande aufge- 
tragen. Fir eine diinne Gipsplatte ist die Ab- 
strahlminderung sehr gro, wahrend fiir eine 


104B ] aad saa 


500 1000 5000 10000Hz 


——= Frequenz 


Abb. 8. Beispiel fiir das-zu_erwartende Abstrahlverhalten 
verschiedener Wandflichen, die zu freien Biege- 
wellen angeregt sind. GréBe der Wiinde zu 3 x 3 m? 
angenommen. 


Kurve a: 12 em Vollziegel, 
Kurve b: 7 cm Porenbeton, 
Kurve c: 1 cm Glas, 
Kurve d: 1 cm Gipsplatte. 
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12-cm-Ziegelwand bereits keine nennenswerte 
Wirkung mehr zu erwarten ist. 

Die Kenntnis dieser Zusammenhange ist fiir 
die praktische Ausbildung von zweischaligen 
Wanden und Decken von Bedeutung. Uber die 
Ergebnisse einer experimentellen Uberpriifung 


wird demnachst gesondert berichtet werden. 
(Eingegangen am 21. November 1952.) 
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BOOKREVIEW 


G. W. Swenson jr., Principles of modern acoustics. 
D. van Nostrand Co. Inc., New York 1953, 222 pages, 
63 fig., 15x23 em’, $ 4.00. 


This textbook is obviously meant as an introduction in 
acoustics for physical and electrotechnical students. It 
covers in a concise form the whole field of acoustics and is 
rather conventionally subdivided into the following chapters: 
I. Oscillations in lumped-parameter systems, II. The vibrat- 
ing string, [1I. The vibrating membrane, [V. Acoustic waves, 
V. Radiators, VI. Acoustic circuits, VII. Architectural acous- 
tics, VIII. Speech and hearing, IX. Approximate methods 
for acoustic problems. Complex calculation and the imped- 
ance concept are used throughout the book. References to 
current literature are frequently made. Every chapter ends 
with several instructive problems. 

Mathematics are used when necessary, but much stress is 
always laid on the physical understanding of the problem in- 
volved. The book is easy to read and yet gives much infor- 
mation. Analogies between mechanical, acoustical and elec- 
trical systems are frequently discussed. Chapter IX (Appro- 
ximate methods) is very useful, on the one hand as a descrip- 
tion of possible mathematical tools for achieving approxi- 


mate solutions, on the other hand to emphasize that rigorous 
solutions of many practical problems cannot be given. 

The MKS-system of units is used throughout, thus break- 
ing with the bad habit of using either a mixture of systems 
or an impractical one. The reviewer congratulates the author 
on this initiative. 

Several shortcomings that need improvement are in our 
opinion: a) p. 13, unit of acceleration m/s/s; p. 14, viscous 
friction in N/m/s (better m/s? and Ns/m resp., or m/s? and N 
per m/s); b) p. 23, eq. (34) f is not (OW/dx)p but (OW/dx)o; 
c) p. 84, the name “‘intrinsic resistance”’ seems bad, why not 
the generally used name “‘wave impedance” as on p. 43; 
d) p. 89, Fig. 4-3, it seems almost certain that an interfero- 
meter with a slot as shown would not function well; p. 152, 
the absorption is not a maximum when 6 = 0; e) p. 136 and 
137, eq. (10) is not a “trefinement” (HELMHOLTZ’ solution is; 
see also CRANDALL, Appendix A), furthermore how can a 
radius be compared with a frequency (p. 137)? f) not Kg M, 
but kg m (p. 137, 213)? 

Notwithstanding these (and possibly some other) faults 
or slips of the pen we think the book to be a very valuable 
introduction in the field. Print, paper and appearance are 
excellent. C. W. Kosten. 


